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ABSTRACT 

A  theoretical  model  is  developed  for  making  computerized  fore- 
casts of  mixed  layer  depth.     An  empirical  equation  relating  significant 
wave  height  to  layer  depth  is  used  to  account  for  mechanical  mixing. 
Heat  exchange  parameters  are  introduced  into  the  model  to  produce 
convective  mixing.      The  u  and  v  components  of  surface  current  are 
used  to  investigate  the  effect  of  convergence  and  divergence.      The 
layer  depths  forecast  by  the  model  are  compared  to  actual  layer  depths 
and  to  the  Fleet  Numerical  Weather  Central  layer  depth  analysis. 
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I.     INTRODUCTION 

With  the  increased  emphasis  on  antisubmarine  warfare  and  the 
development  of  various  sonar  range  prediction  models,    it  has  become 
increasingly  important  that  an  accurate  prediction  of  the  ocean  thermal 
structure  be  raade  available.      Changes  in  ocean  thermal  structure, 
which  occur  mainly  in  the  surface  layer  and  the  thermocline,   are  of 
major  importance  in  deterinining  variations  in  the  vertical  sound 
velocity  profile.      The  sound  velocity  field  determines  the  refraction 
pattern  and  the  existence  of  surface  ducts,    sound  channels  and  other 
sound  transmission  paths.      Therefore,    better  data  on  the  thermal 
structure  will  increase  the  reliability  of  predicting  sonar  perforjnance. 

The  inost  important  thermal  structure  parameter  is  the  vertical 
extent  of  the  surface  duct,    or  the  mixed  layer  depth.     A  surface  duct 
will  exist  if:     (1)  the  temperature  increases  with  depth,    or  (2)  an 
isothermal  layer  extends  downward  from  the  sea  surface.     In  case 
(1),    sound  velocity  increases  with  depth  as  temperature  and  pressure 
increase;  in  case  (2),    there  is  no  temperatixre  gradient  and,    in  the 
absence  of  a  salinity  gradient,    pressure  cavises  sound  velocity  to 
increase  with  depth.      The  greater  the  depth  of  the  mixed  layer,    the 
greater  is  the  difference  between  the  sound  velocity  at  the  surface  and 
at  depth.      The  nuinber  of  sound  rays  which  are  trapped  in  the  layer  is 
therefore  greater.     Also,    the  deeper  the  layer,    the  fewer  the  surface 
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reflections  which  are  required  to  reach  a  given  range,    with  a  resultant 
decrease  in  propagation  losses.     With  increasing  layer  depth,    the 
detection  range  increases  for  a  target  at  its  best  depth  to  avoid 
detection.     Thus  it  is  evident  that  the  mixed  layer  depth,    to  a  great 
extent,    determines  the  sonar  ranges  which  will  prevail.      These 
variations  in  ray  propagation  with  layer  depth  are  illustrated  in 
figure   1. 

The  mixed  layer  depth  forecasting  scheme  presently  in  use  at 
Fleet  Nuinerical  Weather  Central  (FNWC)  is  based  primarily  on 
climatology  and  analysis  of  previously  reported  conditions.     Each  day, 
the  previous  forecast  field  is  adjusted  towards  climatology  and  then  is 
modified  by  turbulent  mixing  due  to  wave  action.      This  method  of 
predicting  mixed  layer  depth  is  rather  elementary:     the  use  of  this 
model  gives  very  little  improvement  over  a  prediction  based  on 
climatology  or  persistence  (forecast  conditions  will  be  the  same  as 
present  conditions).     An  exception  occurs  when  wind  mixing  causes  a 
significant  change  from  climatology;  under  these  conditions  a  fair 
correlation  exists  between  the  forecast  and  actual  conditions.     A  more 
complete  explanation  of  the  method  presently  in  use  at  FNWC  is   con- 
tained in  Appendix  A. 

The  purpose  of  this   research  is  to  produce  a  computerized  mixed 
layer  depth  model  for  FNWC  which  will  make  a  prediction  of  layer 
depth  based  on  environmental  conditions  to  replace  the  present  one 
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Figure  1.      Variation  of  sound  ray  propagation  with  layer  depth. 
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which  is  strongly  dependent  on  climatology.      The  prediction  method 
will  be  considered  improved  if  the  mixed  layer  depth  produced  proves 
to  have  a  smaller  percent  error  when  compared  to  the  verification 
layer  depth  than  the  layer  depth  produced  by  the  FNWC  analysis.     A 
further  check  as  to  the  accuracy  of  the  prediction  can  be  raade  by  com- 
paring the  standard  deviation  of  the  forecast  layer  depth  from  the 
verification  layer  depth  to  the  standard  deviation  of  the  FNWC  method: 
25  feet  in  the  summer  and  80  feet  in  the  winter. 
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II.     PREVIOUS  MIXED  LAYER  DEPTH  STUDIES 

Numerous  studies  relating  to  various  aspects  of  the  mixed  layer 
depth  have  been  conducted  at  the  U.    S.    Naval  Postgraduate  School. 

Clark  [4]  directed  her  attention  to  the  upper  layers  of  the  atmos- 
phere and  such  conditions  that  might  possibly  be  associated  with 
variations  in  the  teinperatvxre  structure  of  the  ocean.     She  presented 
two  hypotheses:     (1)  the  upper  mixed  layer  of  the  ocean  and  the  layer 
of  air  in  iminediate  contact  with  the  ocean  surface  should  be  treated 
as  a  unit  which  is  affected  by  meteorological  parameters  above  this 
unit,    (2)  there  exists  in  the  upper  atinosphere  a  'mirror  image'  level 
where  fluctuations  in  wind  speed  closely  depict  oscillations  in  the  depth 
of  the  mixed  layer.      Clark  reached  the  conclusion  that  a  'mirror  image' 
level  does  exist  in  the  atmosphere  at  an  altitude  of  6000  meters. 

Geary  [lO]  developed  a  model  of  mixing  by  wind  generated  wave 
motion  and  coixipared  it  to  Laevastu's  [18]  and  Neumann's  [23]  methods. 
He  concluded  that  Laevastu's  method  gave  the  closest  approximation  to 
an  upper  limit  on  the  observed  mixing.     He  also  stated  the  need  for  the 
development  of  a  convective  mixing  model. 

Edgren  and  MacPherson  [8]  used  multiple  regression  analysis  on 
various  paraineters   (primarily  meteorological)  to  examine  the  physical 
processes  causing  increases  and  fluctuations  in  the  mixed  layer  depth. 
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It  was  determined  that  convection  is  the  process  which  causes  the 
seasonal  increase  of  the  layer  depth  during  the  cooling  season  and  that 
short  term  fluctuations  of  the  mixed  layer  are  due  primarily  to  internal 
waves. 

Luskin  [21]  developed  a  simple  model  for  raixed  layer  depth  change 
by  thermohaline  convection  and  outlined  a  method  for  adapting  the  con- 
vection model  for  use  as  a  forecast  tool.     He  pointed  out  the  need  to 
develop  a  heat  budget  model  for  forecasting  convective  mixing. 

Lambright  [20]  investigated  the  apparent  random  oscillation  of  the 
mixed  layer  depth  and  by  the  nature  of  its  energy  spectrum  and  its 
correlation  with  tidal  activity,    determined  the  source  of  some  com- 
ponents of  this  oscillation. 

Davis  [5]  investigated  Bulgakov's  model  relating  the  depth  of  con- 
vective mixing  to  density  and  salinity  changes  in  the  upper  layer  of  the 
ocean.     He  sxiggests  the  use  of  this  model  in  taking  into  account  the 
convective  mixing  due  to  increased  density  caused  by  mixing  two  or 
more  layers  of  the  ocean  during  the  heating  season. 

Hancock  [12]  verified  the  Fleet  Numerical  Weather  Central 
analyses  and  forecasts  of  niixed  layer  depth.     The  model  verified  in 
this  study  included  convergence/divergence  computations.      The  con- 
vective mixing  influence  was   calculated  from  a  change  in  sea  surface 
temperature. 

Ciboci  [3]  verified  Klazeika's  [22]  method  for  Ocean  Weather 
Station  "P".     It  was  found  to  be  successful  during  the  heating  season 
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followed  by  a  rapid  decline  as  the  cooling  season  began.     He  also 
determined  that  the  forecasting  curves  developed  by  James  [15]  are 
easier  to  use  than  Mazeika's  and  provide  a  more  accurate  mixed  layer 
depth  forecast. 

Kelley  [16]  studied  a  model  based  on  Kitaigorodsky 's  application 
of  similarity  theory  and  modified  by  McDonnell  to  forecast  mixed 
layer  depth. 

In  this  present  work,    an  attempt  is  made  to  use  various  tech- 
niques previously  developed  to  produce  a  forecasting  method  which 
takes  into  accotint  heat  budget  calculations  in  computing  the  effect  of 
convective  mixing.      The  need  for  a  model  of  this  type  has  been  stated 
many  times  and  this  study  is  a  further  step  in  the  continuing  quest  to 
develop  an  accurate  forecasting  model  for  mixed  layer  depth. 
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III.      FACTORS  AFFECTING  MLD 

The  ocean  may  be  considered  to  consist  of  three  distinct  layers: 

(1)  the  mixed  upper  layer,    influenced  by  atmospheric  conditions; 

(2)  a  zone  of  rapidly  decreasing  temperature  called  the  thermocline 
also  influenced  by  atmospheric  conditions,    and  (3)  a  region  of  slower 
changes  extending  from  the  thermocline  to  the  ocean  bottom. 

It  is  a  feature  of  the  mixed  upper  layer,    specifically  the  depth  at 
which  this  layer  ends  and  the  thermocline  begins,    that  is  of  priinary 
interest  in  this  work.      This  depth  is  termed  the  mixed  layer  depth  (or 
simply  layer  depth)  and  is  defined  as  the  depth  of  maximium  tem- 
perature or  the  extent  of  the  isothermal  layer.      The  term  isothermal 
is  defined  to  include  all  positive  gradients,    and  negative  gradients  up 
to  and  including  0.  3    F/lOO  feet  [15].      This  definition  is   equivalent  to 
the  definition  of  sonic  layer  depth,    the  depth  at  which  maximum  sound 
velocity  occurs  in  the  layer  above  the  therinocline  [2],      Thus  the  two 
terms  may  be  used  interchangeably. 

A  number  of  methods  are  available  for  predicting  the  thermal 
structure  of  the  oceans.     These  involve  calculations  of  net  heat  ex- 
change across  the  air- sea  boundary,    convective  and  wind  mixing, 
continuity,    and  advective  processes.      Each  procedure  may  be  utilized 
separately,    but  the  inost  accurate  forecasts  are  generally  made  by 
considering  the  collective  effects  of  all  processes. 
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A.  NET  HEAT  EXCHANGE 

The  sea- surface  temperature  and  the  thermal  structure  imme- 
diately below  the  surface  and  its  variations  are  greatly  influenced  by 
the  heat  exchange  between  the  sea  and  the  atmiosphere.     It  is  this  energy 
exchange  at  the  air- sea  interface  which  forms  the  basis  for  the  heat 
budget  method  of  forecasting  near-surface  thermal  conditions.      The 
net  heat  gain  or  loss  at  tlie  ocean's  surface  for  the  period  of  the  fore- 
cast is  computed  from  known  or  forecast  values  of  sun's  altitude, 
cloud  cover,    humidity,    sea-surface  temperature,    wind  speed,    relative 
vapor  pressure,    and  albedo.      Of  these  paraineters,  only  the  sun's 
altitude  is  a  truly-known  quantity,    while  the  others  are  forecast  with 
varying  degrees  of  accuracy.     Heat  budget  computations  are  also  of 
use  in  solving  other  problems,    such  as  ice  prediction,    prediction  of 
transient  thermoclines  and  possibly  the  improvement  of  mediuin- range 
weather  forecasts  over  the  sea. 

Once  the  total  heat  exchange  has  been  calculated,    the  gain  or  loss 
of  energy  must  be  distributed  throughout  the  water  column.     Qual- 
itatively,   heat  gains  often  are  associated  with  decreased  layer  depth 
while  losses  accompany  layer  depth  increases. 

B.  MIXING 

The  primary  types  of  mixing  considered  in  the  formation  of  the 
mixed  layer  are:  (1)  instability  (or  convective)  mixing  produced  by 
the  sinking  of  dense  water;  (2)  mixing  as  the  result  of  convergence  or 
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divergence  causing  the  sinking  or  rising  of  water;  and  (3)  mechanical 
(or  forced)  mixing,    a  turbulent  transfer  of  momentum  from  one  level 
to  another  by  the  combined  action  of  wind  waves  and  associated  wind 
currents.     All  three  types  of  mixing  often  occur  simultaneously. 

After  calculating  the  magnitude  of  the  heat  gain  or  loss,    convective 
mixing  and  mechanical  mixing  must  be  considered  to  determine  heat 
distribution.      For  a  net  heat  gain,    implying  a  negative  temperature 
gradient  in  the  surface  layer,    wind  raixing  becomes  the  doininant 
factor  in  computing  the  mixed  layer  depth. 

Although  wind  is  the  direct  force  involved  in  mechanical  mixing, 
mixing  in  the  water  column  is  accomplished  throtigh  the  action  of 
secondary  phenomena-  wind  waves  and  pure  wind  currents.     Swell 
can  also  contribute  to  mechanical  mixing  when  combined  with  the 
effects  of  wind  waves  and  currents.      The  inixing  by  waves  is  due 
primarily  to  the  presence  of  waves  of  different  period  and  height  and 
to  the  turbulence  caused  by  breaking.      In  the  majority  of  cases,    the 
seas  are  not  fully  developed,    and  wind  speed  alone  cannot  be  used 
directly  in  the  determination  of  mechanical  niixing  by  the  sea.     In 
determining  the  thickness  of  a  raixed  layer  created  by  wave  action, 
the  surface  wave  characteristics  are  first  deternained  and  then  a 
depth  is  computed  where  particle  motion  and  the  resulting  mixing  are 
negligible. 

The  mixed  layer  depth  due  to  mechanical  mixing  will  depend  on 
the  wave  model  used.     Since  the  singular  wave  forecasting  model 
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(one  v.'hich  relates  wave  height  to  wind  speed  through  a  simple  exponen- 
tial relationship)  is  presently  in  use,    only  the  simplest  relationship 
between  wave  height  and  mixed  layer  depth  is  justified  [19]. 

C.  CONVECTIVE  MIXING 

Although  mechanical  inixing  is  an  important  parameter  in  deter- 
mining the  mixed  layer  depth,    it  is  not  the  only  one.      In  fact,    it  may 
at  times  be  entirely  inasked  by  other  factors,    such  as  convective 
mixing  and  convergence/divergence. 

When  the  net  heat  exchange  is  from  ocean  to  atmosphere  (negative 
heat  exchange),    convective  mixing  is  the  dominant  factor  in  determining 
the  mixed  layer  depth.      The  upper  layers  of  the  ocean  becoine  unstable 
due  to  cooling  and  sink  until  they  reach  a  stable  level.     This  convection 
will  result  in  isothermal  mixing  of  the  water  to  the  level  of  stability. 

After  determining  the  effects  of  convective  mixing,    the  wind- 
mixed  layer  depth  is  computed.     If  wind  mixing  is  deeper  than  convec- 
tion,   the  column  is  mixqd  isothermally  to  the  wind-mixed  depth.     If 
the  convective  layer  is  deeper,    no  further  modification  is  required. 

D.  CONVERGENCE/DIVERGENCE 

Pure  wind  currents  may  exert  another  influence  on  the  layer  in 
addition  to  the  part  they  play  in  forced  mixing.     In  an  area  in  which 
there  is  a  wind- created  divergence  of  surface  water  (waters  horizon- 
tally advected  apart)  the  layer  depth  decreases  since  cooler  water 
from  intermediate  depths  must  well  upward  to  replace  the  divergent 
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surface  waters.     Similarly,    a  wind-produced  convergence  (waters 
flowing  together)  will  cause  an  increase  in  the  layer  depth  due  to  the 
sinking  of  the  accumulated  surface  water. 

E.  ADVECTION 

Advective  changes  to  the  thermal  structure  are  the  most  difficult 
to  define  with  any  degree  of  accuracy.      The  degree  to  which  advection 
will  affect  the  thermal  structure  is  dependent  upon  the  length  of  the 
forecast  period,    the  local  circulation,    the  horizontal  temperature 
distribution,    and  the  wind  conditions.     In  areas  where  advection  is 
determined  to  be  of  consequence,    primary  consideration  is  given  to 
drift  cxirrents  since  their  effects  tend  to  be  dominant.      Permanent 
currents  are  assumed  to  be  geostrophic  (determined  essentially  by 
the  temperature  distribution  if  there  is  negligible  salinity  contribution) 
and  thus  produce  no  heat  advection. 

F.  INTERNAL  WAVES 

A  primary  factor  which  limits  the  predictability  of  the  instanta- 
neous mixed  layer  depth  and  verification  of  the  forecast  is  the  presence 
of  internal  waves  on  the  interface  between  the  mixed  layer  and  the 
thermocline.      Characteristically,    these  waves  have  large  vertical 
displacements  at  the  bottom  of  the  mixed  layer.     While  similar  to  sur- 
face waves,    internal  waves  lack  the  turbulent  interchange  of  momentum 
(a  result  of  wave  breaking)  possessed  by  surface  waves.      Thus  physical 
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mixing  due  to  internal  waves  is  not  great  unless  these  become  break- 
ing waves.      Therefore,    internal  waves  generally  cause  a  periodic  rise 
and  fall  of  the  mixed  layer  depth  but  do  not  seriously  affect  the  average 
position  of  the  bottom  of  the  mixed  layer.      For  this   reason,    as  well  as 
the  lack  of  sufficient  knowledge  concerning  the  spectral  distribution  of 
internal  waves,    they  were  not  considered  in  this  study. 
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IV.     PROCEDURE 

A.      DATA 

In  the  development  of  this  mixed  layer  depth  forecast  procedure, 
it  was  decided  that  available  bathythermograph  (BT)  reports  would  be 
used  to  obtain  the  actual  mixed  layer  depth  (MLD)  for  use  as  a  verifi- 
cation of  the  predicted  value.      For  this  reason,    the  investigation  is 
concentrated  in  the  vicinity  of  permanent  ocean  station  vessels  on  the 
supposition  that  these  locations  would  have  the  highest  density  of  BT 
data  available. 

It  was  also  desired  to  check  the  forecast  model  during  different 
seasons;  the  study  is  confined  to  data  from  the  months  of  January  and 
September  since  these  are  quasi-stationary  periods  when  mean  physical 
parameters  are  relatively  invariant. 

First,    a  two  degree  la;titude-longitude  square  around  each  of 
eight  ocean  station  vessels  was  checked  for  the  availability  of  BT  data 
for  ten-day  periods  in  January   1972  and  September   1971.     Due  to  data 
limitations,    the  areas  around  Ocean  Station  Vessels   (OSV)   "E",    "N", 
and  "P"  were  selected  as  the  locations  for  use  in  developing  the  pro- 
cedure.    The  periods  used  were  10  to  20  September   1971  and   10  to  20 
January   1972.     Since  these  positions  were  not  on  a  grid  point  for  the 
standard  FNWC  analysis  grid  network,    the  data  fields  were  inter- 
polated to  the  OSV  locations  using  linear  interpolation  between  the 
four  surrounding  grid  points. 
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For  the  periods  selected,  the  FNWC  fields  of  analyzed  sea  surface 
temperature,    forecasted  total  heat  exchange  for  the  next  24  hours, 
"u  and  V  components  of  surface  current,    analyzed  sea  temperatiare  at 
600  feet,    and  forecasted  wave  and  swell  height  were  obtained.      These 
fields  were  the  input  data  to  the  forecast  procedxire.     A  brief  descrip- 
tion of  each  field  is  contained  in  Appendices  B  through  F. 

The  order  followed  in  producing  the  forecast  can  be  seen  in  figure 
2.      The  computer  programs  developed  to  produce  the  MLD  forecast  are 
included  as  Appendix  H. 

The  first  step  in  the  process  is  to  input  an  initial  layer  depth. 
This  depth  will  be  extracted  from  a  BT  taken  24  hours  prior  to  the 
desired  forecast  time.      The  next  step  in  the  initialization  phase  is  to 
input  the  data  to  be  used  in  the  forecast  computation. 

B.      FORCED  MIXING  (HEAT  EXCHANGE  POSITIVE) 

The  actual  forecast  procedure  starts  with  the  determination  of 
whether  the  net  heat  exchange  for  the  next  24  hours  is  positive  or  neg- 
ative.    If  the  net  heat  exchange  is  positive,    the  layer  depth  due  to  forced 
wind  mixing  is  coraputed.      Three  methods  of  determining  the  deptli  of 
turbulent  mixing  are  tested  to  determine  which  will  give  the  best  results. 
The  first  raethod  is  an  empirical  formula  based  entirely  on  the 

significant  wave  height  (H,  /^),    suggested  by  Laevastu  [19]: 

MLD  =  12.5  H^^^  (1) 

The  significant  wave  height  used  in  this  formula  is  actually  the  com- 
bined wave  height  (CH)  due  to  both  wave  height  (WH)  and  swell  height 
(SH),    calculated  as  follows: 
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Figure  2.     Mixed  layer  depth  model  flow  diagram. 
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CH  =YWH^+  SH^  (2) 

The  second  method  of  forced  mixing  investigated  is  the  one 
presently  in  use  at  FNWC.      This  also  is  based  primarily  on  the  com- 
bined wave  height  (CH),    with  a  slight  correction  for  the  magnitude  of 
the  temperature  gradient  between  the  surface  and  600  feet  [9]: 

MLD  -  CH  [   11  -  0.  KSST-   T^qq)  J        ^ST  ^  T^^^ 

(3) 
MLD  =  11-CH  SST  <  T^QQ 

The  final  scheme  tested  is  that  suggested  by  James  [15].      The 

first  step  in  this  method  is  to  calculate  the  absorption  of  radiation  by 

layers.     Absorption  coefficients  by  layers  for  various  water  types 

(clear  oceanic,    average  oceanic,    average  coastal)  are  given  by  James, 

For  the  purposes  of  this  study,    average  oceanic  turbidity  is  assuined. 

This  is  a  valid  assumption  since  this  model  is  not  intended  for  use 

in  coastal  waters.      For  oceanic  water,    with  a  net  heat  exchange  at 

-2 

the  surface  of  500  cal-cm      ,    the  resulting  temperature  difference 

between  clear  oceanic  and  average  oceanic  water  in  the  0-20  foot 

o 
layer  is  0.  08    F.      This  difference  is  not  of  sufficient  magnitude    to 

justify  the  additional  computer  time  required  to  miake  a  determination 

of  water  type.     In  addition,    absorption  will  be  computed  for  only  the 

0-20  foot  layer.      This  is  done  to  decrease  the  number  of  computations 

needed  and  is  considered  valid,    since  an  input  at  the  surface  of 

-2 
500  cal-cm        is  required  to  produce  a  change  in  temperature  of 

o 

0.1    F  in  the  20-40  foot  layer  while  this  same  energy  input  at  the 
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o 
surface  would  cause  a  change  of  1.  3    F  in  the  0-20  foot  layer.     It  is 

the  gradient  in  this   shallow  layer  which  inust  be  overcome  by  inechanical 
mixing.      The  average  change  in  temperature  in  the  layer  due  to  the 
surface  heating  is   calculated,    and  is  then  added  to  the  original  tem- 
perature at  the  ten-foot  level.     A  linear  gradient  is  assumed  through 
this  point  between  the  surface  and  the  original  trace  at  20  feet. 

The  next  step  is  to  compute  a  mixing  parameter  (M    )  which  is  a 

w 

function  of  significant  wave  height  (H^  73)  and  the  period  of  maximum 
energy  (Tj^^x)  U^]: 

Mw  -  T^ax^i/s       •  (4) 

The  T     g^^  in  this  case  is  the  period  of  maximum  energy  which  would 

be  present  for  a  fully  developed  sea  having  the  given  significant  wave 
height.      The  combined  wave  height  (equation  2)  is  again  used  in  this 
case  as  the  significant  wave  height.      Once  this  mixing  parameter  has 
been  computed,    it  is  used  with  the  temperature  gradient  of  the  0-20 
foot  layer  to  compute  the  layer  depth  [15]: 

Ki  -K?    ATM 

LD  ^         (  1  -  e       "^  w)  (5) 

AT 
where  Ki  and  K     are  constants. 

C.      CONVECTIVE  MIXING 

For  the  situation  of  a  net  heat  loss  to  the  atmosphere,    layer 
depth  change  due  to  convective  mixing  must  be  calculated.      This  is 
accomplished  using  a  method  set  forth  by  James  [15].      The  first  step 
is  to  calculate  the  convective  mixing  parameter,    M    : 
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M 


AQi (6) 


c         C,  e        ^  T 
P  \  w 

-7 
where  Q    is  the  net  heat  loss     in  gram  cal.  -cm      ,    C     is  the  specific 
1  P 

o 
heat  of  sea  water  and  is  assumed  to  be  constant  at  0.  935  cal/giTi/    C, 

P       is  the  density  of  sea  water  and  is  assumed  to  be  constant  at 
\  w 

3 
1.024  gm/cm   ,    and     aT  is  the  average  temperature  gradient  from 

the  bottom  of  the  initial  layer  to  6OO  feet  in     F/lOO  feet.      M     is  then 

c 

combined  with  the  initial  layer  depth  (h    )  to  give  the  mixed  layer 
depth  due  to  convection  as  follows: 


The  next  step  is  to  calculate  the  forced  mixed  layer  depth  using 
one  of  the  methods  described  in  the  previous   section.      The  mechanical 
and  convective  mixed  layer  depths  are  then  compared  and  the  deeper 
one  is  taken  to  be  the  actual  layer  depth.     In  various  references,    it 
has  been  suggested  that  when  convective  mixing  is  present,    mechanical 
mixing  can  be  ignored.      This  suggestion  is  also  tested. 

D.      CONVERGENCE/DIVERGENCE 

The  effect  of  convergence  or  divergence  on  the  layer  depth  is 
then  determined.      This  is  accomplished  through  the  use  of 
the  u  and  v  components  of  surface  current,    which  are  combined  to 
give  a  change  in  layer  depth  due  to  convergence  or  divergence  as 
follows  [12]:  3^ 4 
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LD  =  0.864  (U^+U^-U2-U4+V3+V^-V^-V2)    LD/4L 
where  0.  864  is  a  dimensional  constant  to  change  the  variation  in  the 
layer  depth  from  the  metric  system  to  feet,    and  L  is  the  distance 
between  grid  points.      The  decision  to  work  in  feet  rather  than  meters 
was  made  in  order  to  decrease  the  number  of  calculations  (wave 
heights,    swell  heights,    layer  depths  and  temperature  gradients  are 
all  in  feet).      This  value  for  the  change  in  layer  depth  is  then  added  to 
the  layer  depth  previously  calculated  and  the  result  is  the  predicted 
mixed  layer  depth. 

E.  ADVECTION 

The  advection  of  heat  and  the  change  in  the  thermal  structure  due 
to  this  advection  is  then  computed. 

For  the  areas  under  investigation,    it  was  noted  that  due  to  the 

o 

slight  temperature  gradient  between  grid  points  (maximuin.  of  0.  01    F 

per  nautical  mile)  and  the  small  currents  (maximum  transport  on  the 

order  of  50  nautical  miles  per  day)  the  temperature  change  due  to 

o 

advection  would  be  no  more  than  0.  2    F/day  at  any  time. 

Due  to  the  difficulty  in  describing  advection  accurately  and  the  var. 
iable    contribution  it  makes  in  the  areas  under  investigation,    the 
decision  was  made  to  disregard  advection  in  this   study. 

F.  FITTING  THE  FORECAST  TO  DATA 

At  this  point  in  an  operational  procedure,    the  latest  BT  reports 
are  read  into  the  program  and  the  forecast  field  is  adjusted 
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accordingly  with  respect  to  the  reported  layer  depth  and  the  age  of 
the  BT  report.     In  this   study,    this  step  was  omitted  since  the  actual 
BT  reports  were  the  basis  for  determining  the  accuracy  of  the 
forecast. 
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V.     RESULTS 

A.      COMPARISON  OF  MECHANICAL  MIXING  FORMULAS 

In  order  to  determine  the  mechanical  jnixing  method  to  be  used  in 
the  procedure  for  the  prediction  of  mixed  layer  depth,    calculations  of 
forced  mixing  were  made  using  each  of  the  three  proposed  relation- 
ships.    In  rela!  to  the  complete  forecasting  model  (figure  2),    the 
calculations  in  this  section  refer  only  to  the  forced  mixing  step.      The 
heat  exchange  is  given  as  positive.      The  data  input  are  combined  wave 
height  and  the  teinperature  gradient  between  the  surface  and  600  feet. 
All  remaining  steps  in  the  model  are  ignored.      The  computations  were 
made  for  five  specified  periods  when  mechanical  mixing  was  the 
determining  factor  in  the  forecast  layer  depth.      This  is  the  case  when 
the  heat  exchange  is  positive.      The  combined  wave  height  (CH)  v/as 
also  a  factor  in  selecting  the  period  to  be  investigated;  computations 
were  desired  for  a  range  of  wave  heights.      The  layer  depth  from  each 
relationship  and  the  observed  layer  depth  are  compared  in  Table  I. 
The  percent  error  from  the  actual  conditions  (in  parentheses)  and  the 
combined  wave  height  for  the  forecast  period  are  also  shown. 

For  reference,    the  equations  of  the  three  relationships  are  as 
follows: 

1.       Laevastu 

LD  =  12.  5   CH  (9) 
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2.       FNWC 

LD  -  CH  [11-  0.  1(SST-   T^qq)]  SST  ^T^qq  (3) 

LD  =  11-  CH  SST-==T^QQ 


3.       James 

K 
LD 


t —  (1   -   e       ^  w)  (5) 


AT 
A  comparison  of  the  Laevastu  formula  to  the   FNWC  one  shows  that 

the  FNWC  will  always  produce  a  shallower  forecast.      The  FNWC  con- 
stant (11)  is   12%  smaller  than  the  Laevastu  constant  (12.5)  thus  pro- 
ducing a  forecast  layer  depth  due  to  mechanical  mixing  which,    in  the 
absence  of  a  teraperatxire  gradient  between  the  surface  and  600  feet, 
will  be   12%  less  than  the  forecast  produced  by  the  Laevastu  formula; 
as  the  temperature  gradient  between  the  surface  and  600  feet  increases 
in  magnitude,    the  differences  between  the  layer  depths  forecast  by  the 
two  methods  increases. 

In  comparing  the  Laevastu  results  to  those  of  James  it  was  noted 
that  the  percent  difference  between  the  two  layer  depths  varied  with  the 
wave  height.     With  wave  heights  of  the  order  of  five  feet,    the  difference 
between  the  two  models  was  less  than  10%  (absolute  differences  in  the 
range  1-5  feet).     As  the  wave  height  increased  to  10-20  feet,    the  James 
model  produced  layer  depths  approximately   18%  greater.     With  a  large 
wave  height  (26.  6  feet)  the  Laevastu  model  produced  a  layer  over  50% 
deeper  than  the  Jaines  model,    and  18%  shallower  than  the  observed 
layer  depth. 
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TABLE  I.  Comparison  of  layer  depths  in  feet  produced  by  various 

mechanical  mixing  formulas.     Numbers  in  parentheses 

are  percent  error  froin  observed  conditions.      Combined 

wave  height  (CH)  in  feet  and  net  heat  exchange  (Q    )  in 

cal/cm     are  indicated  for  each  case.     Q     is  positive  in 

n        ^ 
all  cases . 


CH 

Q 
n 

LAEVASTU 

FNWC 

JAMES 

OBSERVED 

1. 

4.3 

281 

54  (51.9) 

41  (62.7) 

50  (54.5) 

110 

2. 

4.7 

67 

59  (    1.7) 

45  (25.  0) 

60  (    0     ) 

60 

3. 

10.4 

25 

131  (69.5) 

112  (74.0) 

160  (62. 8) 

430 

4. 

18.5 

257 

231  (43.5) 

200  (51.  1) 

270  (34. 0) 

409 

5. 

26.6 

308 

332  (18.2) 

287  (29.3) 

160  (60. 6) 

406 
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One  aspect  of  the  three  mechanical  mixing  relationships  which 
should  be  noted  is  that  as  the  net  heat  gain  at  the  ocean's  surface 
becomes  larger,    the  FNWC  and  James  models  produce  shallower  layer 
depths  due  to  the  influence  of  the  increased  thermal  gradient  in  the 
upper  layer.      The  Laevastu  formula  contains  no  consideration  of  heat 
gain. 

Only  in  case  2  (Table  I)  was  an  accurate  prediction  of  layer  depth 
under  the  influence  of  mechanical  mixing  produced  (by  the  Laevastu 
and  James  relations).      Cases   1  and  2  occurred  in  Septei-nber  during  a 
period  of  decreasing  surface  heating  and  low  wave  height  when  the 
layer  depths  were  relatively  shallow  (60  to  160  feet).     In  these  cases, 
the  wind  mixing  was  not  sufficient  to  keep  the  surface  layer  mixed 
below  60  feet  and  after  two  days  of  surface  heating,    a  shallow  transient 
layer  was  formed  (the  layer  at  60  feet  was  destroyed  the  following  day 
by  convective  mixing). 

The  Laevastu  relation  exhibited  fair  results  (18%  error)  in  case  5. 
Cases  3  through  5  occurred  in  January  during  a  period  of  increasing 
surface  heating  and  large  combined  wave  height  when  the  layer  depth 
was  fairly  deep  (400  feet).      For  these  cases,    although  the  surface 
heating  was  large,    the  mechanical  mixing  was   of  such  magnitude  that 
a  shallow  surface  layer  could  not  form. 

It  appears  that  there  is  some  relationship  between  mechanical 
mixing,   heat  exchange,    and  the  initial  layer  depth  which  determines 
whether  a  layer  will  be  formed  which  conforms  to  the  mechanical 


34 


mixing  term  or  if  the  mechanical  mixing  will  act  to  naaintain  the 
original  deep  layer  by  preventing  the  formation  of  a  shallow  layer. 
When  the  conditions  are  such  that  the  mechanical  mixing  acts  to 
eliminate  the  shallow  layer,    none  of  the  mechanical  mixing  methods 
investigated  here  produce  accurate  forecasts  for  the  deeper  layer. 

After  a  consideration  of  the  preceding  factors,    it  was  decided 
to  use  the  Laevastu  formula  for  predicting  inechanical  mixed  layer 
depth  in  the  model  tested  in  this  study.     The  Laevastu  formula  was 
chosen  over  the  FTMWC  on  the  basis  of  the  smaller  error  between  the 
forecast  and  observed  layer  depth.      The  James  formula  was  eliminated 
priraarily  due  to  the  difficulty  in  computerizing  it  (computations  in 
Table  I  were  produced  manually).      This  was  not  the  only  considera- 
tion however,    as  the  accuracy  of  the  two  methods   (Laevastu  and 
James)  was  also  compared.     At  wave  heights  up  to  20  feet,    there  was 
little  difference  in  the  two  while  in  the  fifth  case  (wave  height  of 
26.6  feet)  the  Laevast-u  result  was  much  better  (18%  vs.    60%). 

The  difficulty  in  computerizing  the  James  model  arises  from  the 
fact  that  the  constants   (K     and  K    )  in  equation  5  vary  with  the  tem- 
perature gradient.      One  possible  solution  for  this  is  to  use  a  nomo- 
gram provided  by  James  [15]  to  produce  a  grid  of  layer  depths  with 
the  wind  mixing  paranneter  as  one  ordinate  and  the  teraperature 
gradient  as  the  other.      This  could  be  done  but  the  computer  time 
required  for  the  grid  search  in  each  case  would  be  prohibitive. 
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B.      DAILY  FORECASTS 

Using  the  various  procedures  outlined  in  Section  IV,    mixed  layer 
depth  forecasts  were  produced  for  selected  24_hour  periods.     In 
making  the  forecasts  in  this  and  the  following  section,    each  applicable 
step  of  the  model  (figure  2)  was   considered  (the  path  taken  depending 
on  whether  the  heat  exchange  was  positive  or  negative).      For  the 
forecasts  in  this   section,    the  initial  layer  depth  was  determined  from 
the  BT  report  taken  24  hours  prior  to  the  forecast  time.      The  criterion 
for  the  selection  of  a  forecast  period  was  the  availability  of  BT  data 
for  use  as  the  initial  layer  depth  and  as  a  verification  of  the  accuracy 
of  the  forecast.      The  forecasts  are  shown  in  the  computer  output  sec- 
tion (Appendix  I). 

The  purpose  of  this  phase  was  to  check  the  accuracy  of  the  fore- 
cast over  a  24-ho\ir  period  using  an  accurate  layer  depth  as  an  input. 
The  observed  layer  depth,    forecast  layer  depth  and  the   FNWC  analysis 
of  the  layer  depth  at  the  forecast  time  for  each  forecast  are  shown  in 
Tables  II  thru  V.      The  percent  error  between  the  forecast  and  observed 
layer  depth  and  the  FNWC  analysis  and  the  observed  layer  depth  were 
calculated  for  each  case  and  are  displayed  as  the  value  in  paren- 
theses in  each  table.      The  last  line  of  each  table  shows  the  average 
standard  deviation  of  the  forecast  layer  depth  froiTi  the  observed 
layer  depth  and  the  FNWC  analyzed  layer  depth  from  the  observed 
layer  depth. 
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1.  January,    Station  "N" 

The  results  for  this  forecast  period  are  contained  in  Table  II. 

These  forecasts  were  made  for  a  period  of  substantial  negative  heat 

2  2 

exchange  (-400  cal/cm     to  -200  cal/cm   )  with  combined  wave  heights 

of  6-11  feet.     The  procedure  exhibited  improved  accuracy  in  two  of  the 
four  cases.      For  the  periods  when  the  FNWC  analysis  was  more  ac- 
curate,   the  improvement  was  only  about  2%  in  each  case.      The  max- 
imum error  in  the  procedure  was  7.  9%  while  the  analysis  had  two  values 
with  an  error  greater  than  8%.     With  respect  to  predicting  the  average 
conditions  for  the  period,    the  procedure  exhibited  an  error  of  3.6% 
with  an  average  standard  deviation  of  29  feet  as  compared  to  -7.  5% 
and  53  feet  for  the  FNWC  analysis. 

The   FNWC  analysis  was  less  than  the  observed  conditions 
in  every  case  while  the  procedure  forecast  was  higher  in  all  but  one 
case. 

2.  January,    Station  "P" 

Results  are  shown  in  Table  III.      The  first  two  forecasts 

2  2 

were  for  periods  of  small  heat  exchange  (-90  cal/cm     and  25  cal/cin 

respectively)  while  the  remaining  four  were  for  a  period  of  sub- 

2 

stantial  negative  heat  exchange  (-250  cal/cin    ).      The  combined  wave 

height     for  the  period  was   10-13  feet.      For  this  period,    the  FNWC 
analysis  proved  to  be  more  accurate  in  four  of  the  six  cases  and  close 
to  the  model  (less  than  2%  difference  in  percent  error)  in  the  remaining 


37 


TABLE  II.        Comparison  of  observed  MLD,    forecast  MLD,    and 

FNWC  MLD  analysis   (in  feet)  for  Station  "N",    January 
(numbers  in  parentheses  are  percent  error). 


TIME 

DATE 

OBSERVED 

FORECAST 

FNWC 

00 

12 

530 

552(      4.2) 

483(_    8.9) 

12 

12 

570 

561(-    1.9) 

499(-12.5) 

00 

13 

510 

531(      4.  1) 

502(_    1.6) 

12 

13 

530 

572(      7.9) 

497(-    6.2) 

AVERAGE 

535 

554(      3.6) 

495(-    7.5) 

STD. 

DEV 

29 

53 

38 


TABLE  III.      Comparison  of  observed  ML.D,    forecast  MLD,    and 

FNWC  MLD  analysis   (in  feet)  for  Station  "P",    January 
(numbers  in  parentheses  are  percent  error). 


TIME 

] 

DATE 

OBSERVED 

FORECAST 

FNWC 

00 

13 

360 

400(    11.  1) 

391(      8.6) 

*12 

13 

430 

130(    69.  8) 

404(-    6.0) 

00 

19 

390 

440(    12.8) 

389(-    0.2) 

12 

19 

410 

386(_    5.9) 

404(-    1.5) 

00 

20 

380 

403(      6. 1) 

405(      6.6) 

12 

20 

460 

420(_    8.7) 

412(>10.4) 

AVERAGE 

405 

363(_10.4) 

400(_    1.2) 

STD. 

DEV. 

139 

30 

indicates  positive  heat  exchange 
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two  cases.     The  maximum  error  for  the  FNWC  analysis  was   10.4% 
while  that  for  the  procedure  was  -69.  8%.      The  procedure  failed  to 
produce  a  forecast  which  was  within  5%  of  the  observed  layer  depth 
for  any  case. 

In  forecasting  the  average  conditions,    the  procedure  had  an 
error  of  -10.4%  with  an  average  standard  deviation  of  139  feet  while 
the  FNWC  analysis  had  an  error  of  only  -  1.  2%  and  an  average  standard 
deviation  of  30  feet. 

3.       September,    Station  "N" 

Results  are  shown  in  Table  IV.      The  net  heat  exchange  for 

the  first  two  forecasts  is  positive;  the  heat  exchange  for  the  remainder 

2  2 

of  the  period  is  negative,    increasing  from  -83  cal/cm     to  -200  cal/cm 

during  the  period.      The  combined  wave  height  throughout  the  forecast 

period  is  2-6  feet.     In  comparing  the  two  schemes,    the  results  of  OOZ 

on  the   19th  are  discounted  since,    although  the   FNWC  analysis  had  a 

smaller  error,    neither  was  within  500%  of  the  observed  layer  depth 

(this  value  was  included  in  the  calculation  of  the  average  layer  depth 

and  the  average  standard  deviation).      Of  the  remaining   10  forecasts, 

the  FNWC  analysis  proved  to  be  more  accurate  in  six  cases;  however, 

the  best  result  obtained  was  an  error  of  7.  9%. 

The  procedure  showed  an  error  of  15%  in  forecasting  average 

layer  depth  with  an  average  standard  deviation  of  85  feet.      The  FNWC 

analysis  had  a  3.  8%  error  in  predicting  average  conditions  with  an 

average  standard  deviation  of  52  feet. 
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TABLE  IV. 


Comparison  of  observed  MLD,    forecast  MLD,    and 
FNWC  MLD  analysis  (in  feet)  for  Station  "N",  September 
(numbers  in  parentheses  are  percent  error). 


TIME           DATE 

OBSERVED 

FORECAST 

FNWC 

*00 

13 

110 

53(-    51.8) 

132(   20.  0) 

♦  12 

13 

60 

59(-      1.7) 

143(138.0) 

00 

14 

180 

110(-    38.  9) 

141(-21.7) 

12 

14 

160 

60(-    62.5) 

141(      9.2) 

00 

16 

150 

151(        0.7) 

138(-    8.0) 

12 

16 

180 

132(_    26.7) 

l38(-32.  3) 

00 

19 

20 

191(    855.  0)- 

135(575. 0) 

12 

19 

140 

160(      14.3) 

129(-    7.9) 

00 

20 

170 

28(      83.5) 

137(-25.  3) 

12 

20 

170 

141(      17.  1) 

137(-19.  4) 

AVERAGE 

133 

113(-    15.0) 

138(      3.8) 

STD.    DEV. 

85 

52 

indicates  positive  heat  exchange 


41 


I 


4.  January,    Station  "E" 

Results  are  presented  in  Table  V.      The  net  heat  exchange  for 

2 
the  first  two  forecasts  is  approximately  -400  cal/cm    ;   for  the  remain- 

2 
ing  forecasts  the  heat  exchange  is  approximately  -70  cal/cm     with  the 

exception  of  OOZ  on  the   16th  when  there  was  a  positive  heat  exchange 

2 
of  9  cal/cm    .      The  combined  wave  height  during  the  period  is  3-7  feet. 

For  the  first  forecast  layer  depth,    both  methods  show  a  large  error; 

in  the  remaining  five  cases,    the  procedure  presented  in  this  study  is 

more  accurate  four  times. 

In  predicting  average  conditions,    all  forecasts  were  included 

(the  forecast  for  OOZ  on  the   12th  was  not  considered  in  computing  the 

average  standard  deviation).     The  procedure  exhibited  an  error  of 

36.  7%  with  an  average  standard  deviation  of  80  feet  while  the   FNWC 

analysis  had  an  error  of  6%    and  an  average  standard  deviation  of 

240  feet. 

5.  Summary 

Large  deviations  of  the  results  from  observed  conditions 
occur  primarily  when  the  net  heat  exchange  at  the  ocean  surface 
becomes  positive  for  a  short  period  of  time.     In  this  case,    the  layer 
depth  forecast  by  the  model  is  determined  entirely  by  mechanical 
mixing.     When  this  happens,    the  result  is  the  prediction  of  the  forma- 
tion of  a  shallow  layer  (as  the  12Z  forecast  on  the   13th  in  Table  III 
and  the  OOZ  forecast  on  the  20th  in  Table  IV  show)  but  this  may  not 
occur.      To  correct  this  problem,    some  method  miist  be  devised  to 
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TABLE  V.  Comparison  of  observed  MLD,    forecast  MLD,    and 

FNWC  MLD  analysis   (in  feet)  for  Station  "E",    January 
(numbers  in  parentheses  are  percent  error). 


TIME 

DATE 

OBSERVED 

FORECAST 

FNWC 

00 

12 

110 

910(727.0) 

422(284.0) 

12 

12 

580 

661(    14.0) 

434(-25.2) 

00 

15 

620 

770(    24.2) 

531(-14.4) 

12 

15 

650 

679(      4.5) 

544(-16.  3) 

^:-'  00 

16 

100 

47(    53.0) 

493(393.0) 

12 

16 

670 

666(-    0.  5) 

464(-30,7) 

AVERAGE 

455 

622(    36.7) 

481(      5.7) 

STD. 

DEV. 

80 

240 

indicates  positive  heat  exchange 
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determine  whether  this  predicted  shallow  layer  will  actually  form  or 
whether  the  deeper  layer  will  persist.      This  could  perhaps  be  related 
to  some  value  of  wind  mixing  which  will  cause  turbulence  to  such  an 
extent  that  a  shallow  layer  cannot  form.     Another  method  would  be 
to  set  a  limit  (either  absolute  or  a  percentage)  on  the  decrease  in 
layer  depth  which  would  be  perinitted  to  occur.     When  this  limit  was 
exceeded,    the  forecast  mixed  layer  depth  would  be  forecast  by  persist- 
ence (the  forecast  layer  is  the  same  as  the  present  layer  depth). 

Large  errors  between  the   FNWC  analysis  and  the  observed 
conditions  occurred  due  to  the  deviation  of  the  observed  conditions 
from  climatology.      The  variability  of  the  climatological  data  is 
presently  being  evaluated  by  FNWC  and  others. 

Although  internal  waves  were  disregarded  in  this   study,    it  is 
of  value  to  note  their  possible  effect  on  the  accuracy  of  the  forecast. 
Assuming  an  average  amplitude  of  40  feet  for  the  internal  waves 
affecting  these  forecasts,    the  forecasts  can  be  considered  to  be   100% 
correct  if  they  are  within  +40  feet  of  the  average  observed  layer  depth. 
Applying  this  assumed  internal  wave  amplitvide  to  the  forecasts  yields 
an  accuracy  of  90%  in  the  first  two  cases  (Tables  II  and  III).      For  the 
other  two  periods,    the  improvement  is  negligible,    indicating  that 
perhaps  other  factors  are  present  or  that  larger  internal  waves  are 
present. 

Examination  of  the  average  standard  deviations  indicates 
that  although  the  average  standard  deviation  may  be  large,    both  the 
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procedure  and  the  FNWC  analysis  should  be  able  to  predict  the  average 
conditions  with  soine  degree  of  accuracy. 

C.      CONTINUOUS  FORECASTS 

The  purpose  of  this  phase  was  to  see  if  the  procedure  could  take 
an  actual  layer  depth  (from  the  BT  report  taken  24  hours  prior  to  the 
forecast  time)  as  a  starting  point  and  niake  successive  forecasts  for 
a  period  of  time  using  the  forecast  inixed  layer  depth  as  the  starting 
point  for  each  succeeding  prediction.     It  was  hoped  that  after  running 
the  model  for  a  period  of  time  in  this  way,    the  difference  between  the 
resulting  forecast  layer  depth  and  the  observed  condition  would  be 
minimal.      The  forecasts  are  shown  in  the  computer  output  section 
(Appendix  I).  ^ 

The  criterion  used  in  selecting  the  periods  used  was  the  avail- 
ability of  observed  layer  depths  throughout  the  period  for  use  in  inaking 
a  comparison  of  the  accuracy  of  each  forecast.      The  forecast  MLD, 
observed  ML.D,    and  the  percent  error  between  tlie  two  are  shown  in 
Tables  VII  thru  X. 

1.       January,    Station  "N"  (QOZ  forecasts) 

Results  are  shown  in  Table  VI.      The  procedure  worked  well 
(maximum  error  of  8.  4%)  for  the  first  three  days.      On  the  fourth  day, 
there  was  a  positive  heat  exchange,    thus  causing  convective  mixing  to 
cease  and  the  forecast  layer  depth  was  determined  by  mechanical 
mixing.     Since  mechanical  mixing  was  small  on  this  day,    the  layer 


45 


TABLE  VI. 


Comparison  of  the  forecast  MLD  to  the  observed  MLD 
(in  feet)  and  the  percent  error  between  the  two  for 
Station  "N",    January  (OOZ  forecasts). 


DATE 

FORECAST 

OBSERVED 

%  ERROR 

12 

552 

530 

4.2 

13 

553 

510 

8.4 

14 

554 

520 

^.5 

-15 

43 

540 

-92.  0 

16 

82 

440 

-81.4 

17 

115 

540 

-78.7 

^ns 

46 

520 

-91.  1 

19 

77 

500 

-84.6 

20 

87 

590 

-85.2 

AVERAGE 

234 

521 

-55.  1 

indicates  positive  heat  exchange 
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depth  was  forecasted  to  decrease  drastically  (from  520  feet  to  43  feet) 
while  actually  it  increased  by  20  feet.     Since  this  layer  depth  (43  feet) 
was  used  as  the  initial  layer  depth  for  the  next  forecast,    the  large 
error  (80%)  was  propagated  throughout  the  remainder  of  the  period. 
Positive  heat  exchange  again  occurred  on  the   18th  which  caused  the 
layer  again  to  decrease;  this  stopped  the  growth  which  was  occurring 
due  to  convective  mixing. 

Calculations  were  made  to  determine  the  effect  on  the  ac- 
curacy of  the  forecast  if  for  the  forecast  on  the   15th,    persistence  of 
the  layer  depth  on  the  14th  (554  feet)  were  forecast  rather  than  the 
mechanical  mixing  layer  depth  (42  feet).      The  forecast  procedure  is 
then  followed  as  before  until  the  18th  when  again  the  layer  depth  for 
the  previous  day  is  forecast  to  persist.      The  result  of  this  modified 
forecast  is  that  the  procedure  predicts  the  average  layer  depth  with  an 
error  of  6.  3%  and  an  average  standard  deviation  of    53  feet.      These 
values  compare  favorably  with  the  values   coinputed  for  the  three  days 
on  which  good  forecasts  were  obtained  originally  (6.  3%  and  41  feet). 
2.       January,    Station  "P"  (OOZ  forecasts) 

Results  are  sliown  in  Table  VII.      The  results  for  this  period 
were  very  good  althotigh  there  were  three  days  during  the  period  for 
which  no  BT  observations  were  available  for  comparison. 

This  period  was  a  time  of  strong  wind  mixing  (18  to  25  foot 
combined  wave  height)  in  addition  to  being  an  interval  of  substantial 
convective  mixing.      The  large  mechanical  raixing  on  the   14th  and   15th 
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TABLE  VII. 


Comparison  of  forecast  MLD  to  the  observed  MLD 
(in  feet)  and  the  percent  error  between  the  two  for 
Station  "P",    January  (OOZ  forecasts). 


DATE 

FORECAST 

OBSERVED 

13 

370 

390 

♦  14 

233 



^ns 

336 



16 

336 

330 

17 

340 



18 

348 

430 

19 

357 

390 

20 

369 

380 

ERAGE 

336 

384 

%  ERROR 


5.  1 


1 


-19.  1 
-8.5 
-  2.9 
-12.5 


indicates  positive  heat  exchange 
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prevented  a  drastic  decrease  in  the  forecasted  layer  depth  thus  averting 
the  large  errors  which  occurred  at  Station  "N"  as  previously  noted. 

A  stated  desire  of  this  phase  was  to  have  a  minimal  difference 
between  the  final  forecast  layer  depth  and  the  observed  condition. 
This  was  accomplished  for  this  forecast  period  v/ith  a  difference  of 
only  11  feet  (-2.  9%)  on  the  last  day. 

The  average  standard  deviation  for  this  period  was  48  feet 
with  an  error  of  -12.  5%  in  forecasting  the  average  layer  depth. 
3.       September,    Station  "N"  (OOZ  forecasts) 

Results  are  shown  in  Table  VIII.      The  resvilts  were  unsatis- 
factory with  errors  ranging  from  50%  to  85%  between  the  forecasted 
and  observed  layer  depths.      This  was  due  again  to  light  mechanical 
mixing  which  allowed  the  initial  layer  depth  on  the  first  day  to  decrease 
rapidly  (140  feet  to  29  feet)  while  the  layer  acttially  deepened  20  feet. 
This  period  was  dominated  by  positive  heat  exchange  (3  of  4  days)  with 
convective  mixing  on  the  last  day  being  ineffective  in  producing  an 
accurate  forecast  due  to  the  shallow  initial  layer  depth. 

Again,    the  effect  of  forecasting  persistence  for  the  first 
three  days  with  the  calculation  of  the  layer  depth  on  the  fourth  day 
based  on  the  persisted  layer  depth  was  investigated.      The  errors  in 
this  case  between  the  forecast  and  observed  layer  depths  were  in  the 
range,    13%  to  21%.     In  predicting  the  average  layer  depth,    this  mod- 
ified model  had  an  error  of  2%  and  an  av^erage  standard  deviation  of 
33  feet. 
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TABLE  VIII. 


Comparison  of  the  forecast  MLD  to  the  observed  MLD 
(in  feet)  and  the  percent  error  between  the  two  for 
Station  "N",    September  (OOZ  forecasts). 


DATE 

FORECAST 

OBSERVED 

%  ERROR 

*11 

29 

160 

-81.9 

>:a2 

19 

120 

-84.2 

=:=13 

53 

110 

-51.8 

14 

56 

180 

-68.  9 

AVERAGE 

39 

143 

-72.  7 

indicates  positive  heat  exchange 
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4,  September,   Station  "N"  (12Z  forecasts) 

Results  are  shown  in  Table  IX.      This  period  also  produced 
bad  results  with  the  exception  of  the  forecast  for  the   13th  which  had  a 
-1.7%  error.      Errors  on  the  remaining  days  ranged  from  45%  to  70%. 

The  first  three  days  were  controlled  by  the  positive  heat 
exchange  and  thus,    with  no  convective  mixing  occurring,    the  layer 
depth  was  determined  by  mechanical  mixing.     Since  the  mixing  was 
small  (3  to  5  foot  combined  wave  height)  the  surface  heating  produced 
a  shallow  transient  layer  of  60  feet  after  three  days  and  the  only  good 
forecast.     When  convective  iTiixing  occurred  the  following  day,    the 
transient  layer  was  destroyed  but  the  convective  mixed  layer  depth 
model  in  use  is  such  that  a  drastic  increase  in  layer  depth  cannot  be 
forecasted.      Therefore,    although  convective  mixing  was  the  driving 
force  during  the  final  three  days,    the  forecasted  layer  depth  continued 
to  lag  behind  observed  conditions. 

5.  Summary 

Large  deviations  of  the  model  from  observed  conditions 
occurred  mainly  when  an  amount  of  heat  was  added  and  mechanical 
mixing  followed  for  one  or  more  days  causing  the  formation  of  a  shal- 
low mixed  layer  (as  shown  by  figures  3  through  5  of  Appendix  G);  the 
convective  mixing  model  was  unable  to  quickly  eliminate  the  shallow 
thermocline  upon  the  return  to  a  negative  heat  exchange. 

The  one  series  of  good  forecasts  occurred  for  a  period  when 
both  convective  and  mechanical  mixing  were  substantial.     The  large 
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TABLE  IX.  Comparison  of  the  forecast  MLD  to  the  observed  MLD 

(in  feet)  and  the  percent  error  between  the  two  for 
Station  "N",   Septeinber  (12Z  forecasts). 


DATE 

FORECAST 

OBSERVED 

%  ERROR 

^<11 

39 

130 

-70.0 

^nz 

49 

160 

-69.  3 

*13 

59 

60 

-    1.7 

14 

59 

160 

-63.  1 

15 

72 

130 

-44.  6 

16 

74 

180 

-58.  9 

VERAGE 

59 

137 

-56.9 

indicates  positive  heat  exchange 


52 


i 


mechanical  mixing  prevented  the  formation  of  a  shallow  layer  until 
convective  mixing  again  becanae  donainant. 

As  in  the  daily  forecasts  discussed  in  the  previous  section, 
the  need  is  evident  for  some  method  of  predicting  the  formation  of  a 
shallow  layer  during  periods  when  mechanical  mixing  is  dominant. 
Possible  methods  of  accomplishing  this  are  discussed  in  the  summary 
of  the  previous   section. 

D.   CONVERGENCE/DIVERGENCE  CORRECTIONS 

To  determine  the  effectiveness  of  the  convergence/divergence 
corrections,    the  daily  forecasts  were  studied  to  determine  whether  or 
not  the  correction  improved  the  accuracy  of  the  forecasted  layer 
depth.      The  daily  forecasts  were  selected  in  order  to  avoid  errors 
caused  by  poor  input  data. 

In  determining  whether  the  convergence/divergence  corrections 
were  of  any  advantage,    any  correction  less  than  one  foot  was  deter- 
mined to  have  produced  no  effect.     Also,    no  determination  of  the  effect 
of  the  correction  on  the  accuracy  of  the  forecast  was  made  for  those 
cases  in  which  there  was  a  large  discrepancy  between  the  forecast 
and  observed  layer  depths. 

The  first  forecasts  studied  were  for  Station  "N"  in  September  and 
January.  The  maximum  correction  for  this  set  of  data  was  three  feet. 
Of  15  forecasts  which  were  considered  to  exhibit  satisfactory  accuracy. 
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convergence/divergence  had  no  effect  in  nine  cases.     In  the  remaining 
six  cases,    the  correction  improved  the  accuracy  half  the  time. 

The  second  forecast  period  studied  was  Station  "P"  in  January. 
Four  of  the  six  corrections  for  this  period  v/ere  in  the  range,    9  to  1 1 
feet.     Of  the  six  forecasts,    one  correction  produced  no  effect  while 
three  decreased  the  accuracy  and  two  improved  it. 

The  final  period  investigated  was   January  at  Station  "E".      Four 
of  the  six  corrections  in  this   case  were  in  the  range,    11  to  19  feet. 
Applying  the  correction  increased  the  accuracy  in  one  case  while  it 
was  decreased  in  two  others.      There  were  three  forecasts  in  which 
the  correction  had  no  effect. 

It  is  evident  that  although  there  are  locations  where  convergence/ 
divergence  has  very  little  effect  (Station  "N")  there  are  areas  where 
it  must  be  taken  into  account  (Stations   "P"  and  "E").     It  was  also  noted 
that  of  14  occasions  when  the  convergence/divergence  correction  had 
some  effect,    it  improved  the  accuracy  of  the  forecast  only  in  six 
cases. 
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V I .      CONCLUSIONS  AND  R.ECOMMENDATIQNS 

The  results  of  this   study  indicate  that  the  model  which  was 
developed  will  produce  a  slightly  more  accurate  forecast  than  the 
FNWC  formula  when  there  is  substantial  convective  action  and  the 
conditions  vary  from  climatological  conditions   (Table  II).      Under 
these  conditions  the  model  had  a  maximum  error  of  7.  97o  and  a  stan- 
dard deviation  of  29  feet.     This  is  compared  to  a  maximum  error  of 
12.  5%  and  a  standard  deviation  of  53  feet  for  the  FNWC  formula. 

The  greatest  source  of  error  in  the  model  developed  in  this  study 
is  when  convective  inixing  ceases  for  a  short  time.     When  this  happens, 
the  forecast  layer  depth  is  based  on  the  layer  depth  due  to  mechanical 
mixing  and  this  at  times  forecasts  a  decrease  in  the  layer  depth  of  up 
to  400  feet.      This  is  a  serious  flaw  in  the  model  and  work  needs  to  be 
done  in  this  area.     What  is  actually  happening,    is  the  formation  of  a 
new  shallow  layer  (rather  than  the  deep  layer  being  destroyed)  for  a 
short  period  of  time  with  a  rettirn  to  a  deep  layer  as   soon  as   either 
convective  or  mechanical  mixing  becomes  sufficiently  strong  to  destroy 
this  transient  shallow  layer.      Work  must  be  done  in  this  area  to  deter- 
mine under  exactly  what  conditions  a  shallow  layer  will  form  or  be 
prevented  froni  forming.      Two  methods  which  might  be  used  to  improve 
the  accuracy  of  the  forecast  until  such  time  as  improved  methods  are 
available  are  as  follows: 
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1.  Make  a  persistence  forecast  in  the  event  of  a  forecast  of  a 
drastic  change  in  layer  depth. 

2.  Allow  the  layer  depth  to  decrease  only  a  certain  percent  of 
the  difference  between  the  previous  layer  depth  and  the  fore- 
cast layer  depth. 

In  comparing  the  mechanical  mixing  relationships,    the  Laevastu 
method  (used  in  this  study)  was  more  accurate  than  the  FNWC  method. 
In  the  one  good  forecast  of  mechanical  mixing  investigated,    the  FNWC 
relation  had  an  error  of  25%  (compared  to  1.7%  for  the  Laevastu 
relation)  and  the  next  best  result  obtained  was  an  error  of  29.  3% 
(coinpared  to  18.2%  for  Laevastu's  method).     An  area  of  further  study 
concerns  this  Laevastu  method.     It  is  possible  that  the  constant  should 
be  changed,    perhaps  to  a  second  order  term  in  coiTibined  wave  height. 
This  is  suggested  by  the  fact  that  as  the  combined  wave  height  increased 
(from  5  to  27  feet)  the  error  between  the  forecast  and  observed  con- 
ditions decreased  (from  70%  to  18%). 

It  appears  that  convergence/divergence  is  a  factor,    particularly 
for  Stations   "P"  and  "E".     However,    the  method  used  in  this  study  to 
compute  the  correction  for  convergence/divergence  produced  a  correc- 
tion which  improved  the  accuracy  of  the  forecast  less  than  half  the 
time.     Since  FNWC  is  now  in  the  process  of  improving  their  surface 
current  analysis  model,    this  factor  should  be  investigated  at  a  future 
time. 
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The  basic  conclusion  reached  froin  this   study  is  that  the  James 
method  for  convective  mixing  can  be  comptxterized  and  forecasts  can 
be  produced  which  take  the  heat  budget  into  consideration.      The  ac- 
curacy of  this  iTLodel  is  restricted  by  the  accuracy  of  the  input  data, 
particularly  that  of  the  initial  layer  depth.      Given  an  accurate  layer 
depth,    this  model  has  demonstrated  the  ability  to  make  a  continuing 
forecast  using  forecast  layer  depths  as  input  data  (during  periods  of 
large  negative  heat  exchange)  until  such  a  time  as  a  new  BT  report  is 
available  (Table  VIII). 
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APPENDIX  A 
FNWC  POTENTIAL  LAYER  DEPTH  METHOD 

The  first  step  in  the  FNWC  layer  depth  forecast  method  is  the 
construction  of  the  layer  depth  analysis. 

The  initial  step  in  the  analysis  is  to  adjust  the  12 -hour-old 
potential  layer  depth  field  toward  the  interpolated  climatological  field 
of  the  analysis  time.     This  interpolated  field  is  a  linear  interpolation 
over  two  months;  the  present  and  preceding  month  if  the  day  is  less 
than  15,    or  the  present  and  following  iTionth  if  the  day  is  more  than  15, 
If  the  analysis  is  for  the   15th,    only  the  climatology  of  the  present  m 
month  is  used.      This  adjusted  cliniatological  field  is  the  first  guess 
field. 

Next,    the  mixed  layer  depth  possible  due  to  wave  mixing  alone  is 
computed  from  the  combined  wave  height  analysis  as  follows: 

MLD  =  CH[ll-0.  I(SST-T^qq)]  SST  ^  T^qq 

=  11     CH  ,  SST-e  T, 

600 

The  effect  of  the  term,    0.  I(SST-T^qq),    is  to  reduce  the  depth  of 
mixing  when  a  strong  therinocline  gradient  exists.  ' 

The  mixed  layer  depth  field  due  to  wave  action  is  then  smoothed 
with  a  special  smoothing  function  designed  to  decrease  the  highest 
values  in  the  field  without  increasing  the  values  at  any  point. 
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Next,    the  MLD  field  due  to  wave  mixing  is  compared  to  the 
initial-guess  field  produced  from  climatology  and  the  deeper  of  the 
two  is  chosen  as  the  second-guess  MLD  field. 

The  final  step  in  the  analysis  phase  is  to  adjtist  this   second  guess 
field  to  BT  reports  taken  within  the  past  72  hours.      The  reports  are 
weighted,   based  on  the  age  of  the  observations  at  analysis  time  with 
the  latest  reports  receiving  the  greatest  weight.      One  final  control 
placed  on  the  analysis  field  is  to  force  the  layer  depth  to  be  less  than 
or  equal  to  900  feet  at  all  grid  points.      The  result  is  the  FNWC  poten- 
tial mixed  layer  depth  analysis  field. 

The  forecast  layer  depth  field  is  produced  in  the  saine  manner  as 
the  analyzed  field.      The  analyzed  field  is  adjusted  toward  the  clima- 
tology of  the  forecast  time  using  the  method  described  for  producing 
the  analysis.     The  jmcchanical  inixing  is  then  computed  using  the  24- 
hour  forecast  combined  wave  field.      The  same  steps  of  comparison  and 
smoothing  are  taken  between  the  wind  mixed  field  and  the  clima- 
tological  field.      The  result  is  the  forecast  potential  mixed  layer  depth 
field. 

The  one  difference  in  procedure  between  the  analysis  and  fore- 
cast phase  is  that  the  forecast  field  is  not  adjvisted  to  BT  reports. 

A  severe  constraint  on  this  method  is  the  lack  of  data  which 
results  in  relying  essentially  on  climatology  in  areas  of  sparse  data. 
An  additional  problem  area  is  that  this  forecast  scheme  is  also  subject 
to  errors  in  the  wave  height  forecast. 
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APPENDIX  B 
FNWC  SEA  SURFACE  TEMPERATURE  ANALYSIS 

The  sea  surface  temperature  analysis  is  accomplished  through  a 
scheme  called  Fields  By  Information  Blending  (FIB).     The  analysis 
is  done  on  a  125x125  grid  with  an  output  on  a  63x63  grid. 

Inputs  to  the  analysis  are  sea  surface  temperature  reports  within 
six  hours  of  the  analysis  time  and  surface  temperatures  frora  BT 
reports   received  in  the  past  12  hours.     When  it  becomes  readily 
available,    satellite  data  also  will  be  incorporated  in  the  model. 

The  first  step  in  the  analysis  is  to  compute  the  current  clima- 
tology.     This  climatology  is  then  assigned  a  weight  based  on  the 
gradient  of  temperature  between  grid  points   (highest  gradients  are 
given  the  lowest  weights).     Next,    the  first  guess  field  is  produced  by 
blending  the  previous  analysis  and  the  current  climatology. 

The  next  step  is  to  read  in  and  screen  the  sea  surface  tem- 
perature reports.     A  difference  field  is  computed  by  subtracting  the 
first  guess  field  from  the  field  produced  from  reports.      Each  report 
is  assigned  a  weight  based  on  its  age. 

The  final  sea  surface  temperature  analysis  is  produced  by  blend- 
ing this  difference  field  with  the  guess  field  using  the  assigned  weights, 
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The  basic  drawback  to  this  procedure  is  the  fact  that  1100  to  1200 
reports  are  input  to  produce  an  analysis  over  3969  grid  points.      Ob- 
viously,   some  accuracy  is  lost  through  sraoothing.      The  procedure 
should  become  better  once  satellite  data  are  routinely  available  for 
use  as  input  parameters. 
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APPENDIX  C 
FNWC  SUBSURFACE  THERMAL  STRUCTURE  ANALYSIS 

The  subsurface  temperature  analysis  is  used  for  the  following 
standard  levels:     surface,    100  feet,    200  feet,    300  feet,    400  feet,    600 
feet,    800  feet,    and  1200  feet. 

The  inputs  to  the  model  are:      the    12 -hour-old  teinperature  analysis 
at  standard  levels,    climatology  at  standard  levels,    sea  surface  ten:i- 
perature  analysis,    mixed  layer  depth  analysis,    and  BT  observations 
for  the  past  72  hours. 

The  first  step  is  to  prepare  guess  fields  for  each  standard  level. 
This  is  done  in  two  steps.      First  the  12-hour-old  analysis  is  adjusted 
toward  climatology;  then  this  field  is  adjusted  to  the  sea  surface  tenn- 
perature  and  mixed  layer  depth  analyses  to  produce  the  final  guess 
field. 

The  second  step  is  to  prepare  data  lists  from  reported  observa- 
tions.    Each  report  is  assigned  a  weight  factor  on  the  basis  of  age  and 
then  interpolated  to  find  temperatures  at  all  standard  levels  from  the 
surface  to  the  bottom  of  the  report. 

The  final  step  is  to  adjust  the  fields  to  the  data.      First  the  final 
guess  field  is  adjusted  to  the  data  and  a  light  smoothing  function  is 
applied  to  the  vertical  gradient  between  the  level  in  question  and  the 
level  above.      This   results  in  a  first  analysis  field.     Next,    the  first 
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analysis  field  is  adjusted  to  the  data  and  a  vertical  control  is  applied 
requiring  that  the  temperature  at  the  level  in  question  is  less  than  or 
equal  to  the  temperature  of  the  level  above.      The  final  analysis  for  a 

level  is  produced  after  the  temperature  at  any  level  is   restricted  to 

o  o 

the  range  -2    C  to  35    C. 

The  primary  liinitation  in  this  method  is  the  lack  of  data.     In 

addition,    areas  of  strong  currents  and  areas  of  strong  winds  which 

are  not  forecast  will  produce  errors  in  the  forecast. 
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APPENDIX  D 

FNWC  HEAT  EXCHANGE  MODEL 

The  total  heat  exchange  across  the  air-sea  interface  is  a  com- 
bination of  five  components:    hea.t  absorbed  fronri  the  sun,    solar  heat 
reflected  from  the  sea  surface,    longwave     radiation  from  the  sea  sur- 
face  less  the  incoming  longwave  radiation  from  the  atmosphere, 
sensible  heat  exchange  across  the  air -sea  interface,    and  latent  heat 
exchange  thru  evaporation  or  condensation.      Each  component  was 
calculated  from  a  separate  equation  and  the  results  were  combined  to 
obtain  the  net  heat  exchange.      Following  are  the  equations  used  in  the 
computation  of  each  component: 

1.  Solar  insolation,    Q 

s 

3 

Q     =0.014At^{1.0-0.  0006  C    ) 
s  n  d 

2.  Reflected  solar  insolation,    Q 

r 

Q     =  0.  15Q     -  (0.  OIQ    )^ 
r  s  s' 

3.  Effective  back  radiation,    Q 

b 

0^  =  0^^(1.0-0.765  0)  where, 

Q^b  =  O'T^  [A+4(  '^w  -  '^a  f]        ^,^^^^^^ 

^^       -4 
A  =  0.  261  exp  [-7.  77  x  10        (273-T    )] 

a 

4.  Latent  heat  transfer,    O 

e 

O     =  (0.26+0.  077V)  (0.  98e     -e.)L,     [e     -e^        positive] 
e  '^wa't'-wa-'^  -• 

Og  -  0.  077V{0.  98e^-e^)L^_  t^W^a        negative] 
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5.       Sensible  heat  transfer,    A 

Q,    =  39(0.26  +  0.  077)  (T    _T    )  [T    -T  positive] 

Q^   -  3V(T,„  _   T    )  [T     -T  negative] 

These  five  parameters  are  then  con:ibined  to  give  total  heat  exchange, 

Q   ,   as  follows: 
n 

n  s       r       b       e       h 

These  formulas  were  selected  through  a  coinprehensive  literature 

review  and  were  those  which  promised  the  greatest  reliability  because 

of  performance  in  studies  supported  by  field  testing. 

The  accuracy  of  these  forrnulas  in  synoptic  use  can  be  determined 
directly  only  in  rare  instances  when  a  research  vessel  is  available  to 
make  radiation  measurements.     However,    it  is  possible  to  make  an 
indirect  verification  by  comparing  the  heat  exchange  analyses  with 
actual  measured  changes  which  occur  in  the  ocean  and  atmosphere. 
This  type  of  indirect  verification  has  been  used  for  the  past  several 
years,    and  there  has  been  no  evidence  that  the  FNWC  heat  exchange 
formulas  are  inaccurate  or  in  need  of  change. 

The  errors  which  do  occur  in  the  heat  exchange  computations  are 
probably  caused  by  inaccuracies  in  the  input  analyses  and  forecasts. 
In  general,    these  input  parameters  are  reasonably  accurate  but  errors 
in  meteorological  observations  and  in  the  manipulating  and  averaging 
of  sparse  data  over  large  sea  areas  can  cause  problems.      Cloud 
analyses  have  been  one  of  the  greatest  error  sources.      The  wind  field 
analyses  are  generally  good,    but  a  small  error  in  the  wind  field  at 
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high  wind  speeds  will  cause  large  errors  in  some  heat  exchange 
components,    particularly  in  latent  heat  exchange. 
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APPENDIX  E 
FNWC  SURFACE  CURRENT  ANALYSIS 

Ocean  currents  are  caused  by  the  combined  action  of  several 
physical  environmental  forces.      Thermohaline  circulation,    wind  and 
mass  transport  by  waves,    and  tidal  effects  are  the  most  important 
of  these  forces.     Other  forces  such  as  changes  in  atmospheric  pressure, 
inertia  currents,    current  boundaries,    and  the  coriolis   effect  were 
considered,   but  were  discarded  for  lack  of  significance  or  for  lack  of 
any  forecasting  method.      Tidal  currents  are  important  in  coastal 
waters,    and  have  some  magnitude  in  deep  water.     However,    it  was 
decided  that  for  a  24-hour  prediction,    the  tidal  currents  could  be 
neglected.     The  elimination  process  finally  led  to  the  assumption  for 
the  model  that  the  total  current  flow  could  be  described  by  the  vector 
suin  of  permanent  thermohaline  flow  and  local  wind-driven  flow, 
including  mass  transport  by  waves. 

The  thermohaline  or  permanent  circulation  evolves  from  the 
distribution  of  ocean  temperature  and  salinity.     The  thermohaline- 
type  currents  have  been  studied  for  soiiie  70  years,    and  have  been 
modeled  by  geostrophic  methods.     A  drawback  in  applying  this  geo- 
strophic  method  is  that  it  requires  accurate  information  on  the  dis- 
tribution of  temperature  and  salinity  in  the  water  column  for  which 
little  real-time  salinity  data  is  available.      Therefore,    the  direct 
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application  of  the  geostrophic  method  at  FNWC  has  been  modified  to 
calculation  of  a  permanent  flow  component  based  on  the  local  tem- 
perature structure  alone. 

The  permanent  component  calculation  starts  by  computing  a 
weighted  average  teraperature  from  the  surface  to  the  600-foot  depth. 
Values  for  the  two  levels  are  taken  from  FNWC  synoptic  analyses. 
Using  this  average  temperature,    the  thermal  wind  equation  is  applied 
to  obtain  u  and  v  components  of  tlie  surface  current  caused  by  distribu- 
tion  of  sea  water  temperatures. 

The  second  component  in  the  surface  current  model  is  the  local 
wind-driven  flow.      Coinplicated  calculation  methods  are  available, 
but  they  are  not  well  adapted  to  synoptic  analysis  so  an  empirical 
approach  has  been  used.     A  simple  foriiiula  (a  constant  times  the  wind 
speed)  was  developed  relating  wind  to  surface  currents.     The  constant 
is  adjusted  to  account  for  mass  transport  by  waves.      The  wind  is 
computed  froin  the  FNWC  sea  level  pressure  analysis  and  the  24-hour 
sea  level  pressure  forecast.      The  geostrophic  wind  is  corrected  for 
stability  and  curvature  to  obtain  a  "marine"  wind. 

Thus,    the  FNWC  total  ocean  surface  current  is  the  vector  sum  of 
the  thermal  and  wind-driven  components. 

The  surface  current  model  is  designed  to  produce  an  accurate 
picture  of  current  patterns  over  deep  water.  It  will  do  this,  but  it 
makes  no  pretense  of  accuracy  in  other  areas  and  for  all  possible 
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situations.     The  grid  size  of  about  200  miles  precludes  accuracy  in 
areas  of  marked  therinal  gradients,    such  as  along  the  edge  of  strong 
currents  like  the  Gulf  Stream.     Very  narrow  currents  can  pass  between 
grid  points  without  showing  in  the  surface  ciirrent  analysis.     The  model 
does  not  perforra  well  where  for  those  few  currents  which  are  dom- 
inated by  a  haline  distribution,    as  in  areas  where  there  is  much  melt- 
ing of  ice  or  alternate  periods  of  melting  and  freezing.     A  final  area 
of  model  weakness  is  in  coastal  areas  or  near  island  chains. 
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APPENDIX  F 
FNWC  SINGULAR  ADVECTIVE  WAVE  AND  SV/ELL  MODEL 

The  total  "sea  state"  may  be  described  as  a  spectrum  of  sinusoidal 
waves  of  varying  frequency  and  direction  of  movement. 

The  singular  model  typically  relates  wave  height  to  wind  speed 
through  a  simple  exponential  relationship.     It  is  significant  to  note 
that  singular  models  usually  do  not  maintain  wave  continuity  from  one 
computational  cycle  to  the  next,    nor  is  there  provision  for  energy 
sinks  through  dissipation  due  to  opposing  winds. 

Spectral  models  treat  the  growth,    propagation  and  dissipation  of 
energy  for  a  number  of  frequencies  traveling  in  a  number  of  directions. 
These  models  divide  the  frequency  and  direction  spectrura  into  a  finite 
set  of  bands,    calculate  the  changes  in  energy  for  each  frequency- 
direction  coinponent  and  integrate  over-all  components  to  obtain  energy 
moments  which  statistically  describe  the  total  sea  state.     More 
sophisticated  spectral  models  contain  wave/wave  interaction  terms  in 
addition  to  those  relating  wind/wave  interaction. 

Since  180  components  are  treated  at  each  grid  point,    spectral 
models  are  very  demanding  on  computer  time.      The  singular  advective 
model  presently  in  use  represents  a  compromise  between  computer 
limitations  and  model  completeness.     It  utilizes  the  propagation  of 
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energy  concept  but  only  the  frequency  of  maximum  energy  is  propogated 
in  the  direction  of  the  wind.      (A  single  component  is  selected  as  being 
representative  of  the  total  spectrum.  )     Growth  and  dissipation  are 
achieved  through  comparison  of  advective  heights  with  fully  developed 
heights  from  a  new  wind  field  at  a  later  tiine. 

The  starting  points  for  each  analysis  cycle  are  the  six-hour  old 
wind  wave  and  swell  fields.      The  previous  wind  wave  analysis  is 
advected  forward  in  two  three-hourly  time  steps  with  dissipation  and/ 
or  growth  through  comparison  with  the  current  wind  analysis.      The 
product  of  these  two  advective  steps  is  a  preliminary  calculation  of  the 
new  wind  wave  analysis. 

A  final  wind  wave  analysis  is  obtained  through  the  introduction  of 
synoptic  ship  reports  into  the  preliminary  calculation  of  wave  heights. 

The  forecast  portion  of  the  model  is   simply  a  continuation  of  the 
advective  propagation  outlined  above  except  that  forecast  winds  are 
used  in  the  dissipation  and  growth  of  the  wind  waves. 

The  six-hour-old  swell  trains  are  similarly  advected  forward, 
but  without  consideration  of  any  interaction  with  cxirrent  winds. 
(This  is  a  weakness  which  is  not  inherent  to  true  spectral  models 
where  wind/wave  and  wave/wave  interactions  occur  continually  for 
all  frequencies.  ) 

The  speed  of  propagation  is  proportional  to  swell  period.     Since 
this  model  makes  no  provision  for  the  dissipation  of  swell  energy 


71 


through  swell/wind  interaction,    swell  heights  could  propagate  unchanged, 
To  insure  that  there  will  always  be  some  decay  in  swell  heights  with 
time,    a  decay  coefficient  has  been  added  to  the  advective  equations. 

There  are  three  types  of  limitations  which  contribute  to  inac- 
curacies in  sea/swell  analyses  and  forecasts: 

a.  Wind  errors. 

b.  Wave  observation  errors. 

c.  Modeling  limitations. 

Even  if  model  formulations  were  perfect,    wind  errors  would  lead  to 
incorrect  wave  computations.     Determination  of  the  wind  at  the  surface 
of  the  water  which  actually  controls  wave  growth  is  particularly 
difficult.      The  forecast  winds  used  in  the  forecast  cycle  also  contain 
any  errors  resulting  from  deficiencies  in  the  model  used  to  produce 
the  winds. 

Wave  observation  errors  also  contribute  to  errors  in  the  final 
wind  wave  analysis.     This  is  a  significant  factor  since  almost  all 
wave  height  reports  represent  an  "eyeball"  estimate. 

The  inability  of  any  singular  model  to  define  the  complete  wave 
energy  spectrum  and  all  its  interactions  is  a  major  modeling  limitation. 
The  advective  singular  model  represents  a  significant  iinprovement 
over  earlier  singular  approaches  for  it  insures  wind  wave  continuity 
between  runs  and  treats  energy  dissipation  due  to  opposing  winds. 
However,    wave/wave  interactions  which  influence  both  growth  and 
dissipation  are  not  considered  at  all. 
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APPENDIX  G 


SELECTED  BT  TRACES 


Figures  3   thru  5  show  another  example  of  the  forination  of  a 
shallow  surface  layer.     The  difference  between  this  case  and  the 
preceding  one  is  that  in  this  example  the  layer  is  forraed  not  due  to 
heating  but  due  to  wind  mixing  during  a  period  of  slight  negative  heat 
exchange.     Figure  5  shows  the  return  to  a  deep  layer  following  ^ 
period  of  substantial  negative  heat  exchange. 
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Figure  3.     BT  trace  for  OOZ   18  September,    Station  "N' 
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Figure  4.     BT  trace  for  OOZ   I9  September,    Station  "N". 
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Figure  5.     BT  trace  for  OOZ  20  September,   Station  "N' 
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APPENDIX  H 


PROGRAM    FOR    DAILY    MIXED    LAYER    DEPTH    FORECASTS 


50 
900 


500 
500 


55CC 


50 
INTE 


IIMTE 


INTE 


INTE 


INT 


INTE 


REDU 


DIMENSI 
ISH{^)t 

N=0 

WRITEC  6 
FORMAT { 
N  =  N+1 
IF(N.EQ 
I  F  (  N  .  E  Q 
GC  TO  5 
WRITE (6 
FCRMAT{ 
(CCNT) ' 
READ(5t 
READ(5 , 
READ(5, 
READ(5, 
READ(5, 
READ( 5, 
READ(5t 
READ(5, 
READ(5t 
READ{5, 

1  FORMAT 

2  FORMAT 
0  FCRMAT( 
RPOLATICN 

X  = 
Y  = 

RPOLATE  T 
TL=IT(  L 
TU=IT(3 
TF=TU-Y 

RPOLATE  T 
TSL=ITS 
TSb=ITS 
T5F=TSU 

RPOLATE    H 

.  QKL=IQN 
CNU=IQN 
QNF=QNU 

RPOLATE  W 
WHL=IWH 
WhU=IWH 
WFF=WHU 

RPOLATE  S 
ShL=ISH 
SHU=ISH 
SFF=SHU 

CE  DATA  b 
TC=TF/6 
TSC=TSF 

.  T=((9 
TS=((9. 
QN=ONF/ 
Wh=WHF/ 
SH=SHF/ 
DO  100 
U  (  J  )  =  UJ 
V( J)=IV 


ON  IT(4)  ,ITS(4)  ,IQN(4) ,IU(4) , IV(4)  ,  1WH(4), 
U  {  4  )  t  V  ( ^1  ) 


,501) 

' L S /////// /»27X, 'FORECAST 


MLD  FOR  STATION  PAPA') 


.9) 


GO  TO 
GO  TO 


9999 
5000 


lOL  ) 


.6) 
500 
,50( 

• L* ,////////, 27X, 'FORECAST  MLD  FOP  STATION  PAPA 
) 

50  0 
1  ) 
U 
L) 
1) 
1) 
1  ) 
I) 
2) 
2) 
(41 
(F7 
214 
FA( 


)  KT.KD 
(IT(  1)  r 
(ITS{ I) 
( I QN (  I) 
(1U(  I)  , 
{IV(  I)  , 
{  IWH(  I  ) 
(ISH( I ) 
PL  I 
RLD 
6) 

2) 
) 
XTORS 


1=1,4) 

,1=1,4) 

,  1=1,4) 

1=1,4) 

1=1,4) 

t  1=1,4) 

,1=1,4) 


X    C    Y 


EH 
)- 
)- 
-{ 
EM 
(1 
(3 
-Y 
EA 
(  1 
(3 
-Y 
AV 
(1 
(3 
-Y 
WE 
(  1 
(3 
-Y 
Y 

4. 

/6 

^^TC 

16 
51 
51 
J  = 
(J 
(J 


PE 
X^ 
X- 
TU 
PE 
)- 
J-. 

T 

)- 

)- 

^M 

E 

)- 

)- 

-( 

LL 

)- 

)- 
-( 

SC 


RATU 
(  IT( 
(IT( 
-TL) 
RATU 

x-=(i 

X'^=(I 
TSU- 
EXCH 
X^=(I 
X  =  =  (  I 
UlU- 
hhlG 
X^-'  (I 
X^^  (  I 
WHU- 
HEI 

x-'^ci 

X-^CI 

SHU- 
ALIN 


RE 

1)-IT(2)  ) 
3 )-IT(4)) 


RE     600 


TS(1 

TS(3 

TSL) 

A  NG  E 

QN  (  1 

QN{3 

QNL) 

HT 

WH(1 

WH(3 

i ;  h  L  ) 

GHT 

SH(  1 

SH(3 

SHL) 

G    FACTOR 


-ITS(2)  ) 
-ITS(4) ) 


-IQrM2) 
-ICN(4) 


-IWH(2) 
-I WH(4) 


ISH( 2) 
ISH(4) 


)+l60. )/5. 
SC)+160.)/5. 

• 

2  . 

2. 

1,4 

)/16. 

)/16. 
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100  CONTINUE 
CP=0.935 
RW=  1.024 
DCCN=0  .0 

IFCQN.GE.O.O)    GO    TG    1200 
CGNVECTIVE    MIXING 

DT=(100.--IT-TS)  ) /(PL  I -600.) 
CM=(2.-^QN)/{CP^'^RW-''DT) 
CCCN=(PLI^-2  +  CM)'^"0.5 
MECHANICAL    MI  XING 
1200    Ch=(WH=^--2  +  SH^'=^2)*-0.5 
DM  =  12.5'-^CH 

IF(DM.GE.DCGN)    GO    TO    1300 
DL=DCON 
GC    TO    1000 
1300    DL=DM 
CONV/DIV    CORRECTION 

1000    A  =  U(1 )+U(3 )-U(2)-U(4)+V(3)+V(4)-V(l)-V( 2  J 
GRID    SPACING  'G' 

G= 

B=(A'-''0.864)/(4,=^G) 
DD=B^DL 
PRLD=DL+DD 
WRITE(6,502)    KT»KD 
3000)    DCON 
3  100)    DM 
3200)    DL 
3300)    DD 
3400)    PRLD 
3500)     RLO 

,////, 13Xt'F0RECAST    for    •,I2t'Z 
,23X, • CONVECTIVE    MLD     IS     •tF8.2) 
FORCED    MLD     IS     •  tF8.2) 
LAYER    DEPTH     IS     • ,F8.2) 
CONV/DIV    CGRR    IS     SF8.2) 
FORECAST    ,MLD    IS     •  ,F8.2) 


502 
3000 
31CC 
3200 
3300 
3400 
3500 

9999 


WRITE( 6 

WRITE(6 

WRITE(6 

WRITE(6 

WRITE(6 

WRITE (6 

FORMAT ( 

FORMAT ( 

FCRMAT( 

FORMAT ( 

FCRMAT( 

FORMAT( 

FORMAT ( 

GC  TO  9  000 

STOP 

END 


0 


,23X, 
t23X, 
»23X, 
t23X, 
f23X, 


,12,'  JAN.') 


ACTUAL  MLD  IS  • ,F8.2) 
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PRCGRAM  FOR  CONTINUOUS  MIXED  LAYER  DEPTH  FORECASTS 


501 


9000 


50CC 
5001 

55GC 


DI^'ENSICN 
ISH(4) , U( 
N  =  0 

lvPITE(6,5 
FORMAT{ »L 
STATION  P 
RtAD{5,2) 
PRLD=PL  I 
K  =  N  +  1 
IF(N.EQ.9 
IF(N.EQ.6 
GC  TO  550 
KRITE( 6,5 


IT(4)  ,ITS(4) ,  IQN(4)t  IU(4), IV{4) , IWH(4), 
4) tV(4) 

01) 

',////////. 22X»«C0NTINUCtS    MLD    FCPECAST    FOR 
APA»  ) 
PLI 


TO 
TO 


9999 
5000 


FORKAT( 
(CCNT) • 
READ( 5, 


I 
) 

50 
1) 
1) 


)     GO 
)    GO 

0 

ocu 

',////// //»27X,« FORECAST 


KiLD    FOR    STATION    PAPA 


READ(5, 
PEAD(5, 
RE/iC(5,  n 
PE/^D{5,  1  ) 
READ(5,1) 
READ(5,  L) 
RE/>C(5t1) 
PLI=PRLD 
READ(5, 2) 

1  FCRf/AT     (4 

2  FORMAT     (F 
500    FCRKAT(2I 

INTEPPGIATI  CN  F 
X  = 
Y  = 

INTERPOLATE  TEN! 
TL=IT(  U- 
TU=IT(3  )- 
TF  =  TU-Y^^( 

INTERPOLATE  T EM 
TSL=ITS (I 
TSL=ITS(3 
TSF=TSU-Y 

INTERPOLATE  HEA 

QNL=IQN( 1 

-  GNU=IQN(3 

£NF=gNU-Y 

INTERPOLATE  WAV 
WhL=lWH( L 
lvHtj  =  I  WHO 
UFF=WHU-Y 

INTERPOLATE  SWE 
SHL=ISH(1 
SHU=ISH(3 
SFF=SHU-Y 

REDLCE  CATA  BY 
TC=TF/64. 
T5C=TSF/6 
T=(  (9.'.aC 
TS=((9.^T 
CN=CNF/16 
WH=WHF/ 51 
SH=SHF/51 
DC  100  J= 
U(J)=IU(J 
V(  J)=IV { J 
ICC    CONTINUE 


C)     KT,KD 
(IT(  I) 
(ITS( I 

(  IQN(  I 
(IU(  I) 
(IV(I) 
{ IWH(  I 
( ISH(  I 

RLD 
16) 
7.2) 
4) 

ACTORS 


=  1,4) 
1=1 ,4) 
1=^1,4) 
=  1,4) 
=  1,4) 
1=1,4) 
1=1,4) 


X    &    Y 


PERATURE 
X-M  IT{  1)~IT(2)  ) 
X*(  IT(3  )-IT(4)  ) 
TU-TL) 

PERATURE     600 
)-X*(ITS (1)-ITS(2) ) 
) ~X^(ITS(3)-ITS(4) ) 
===(TSU-TSL  ) 
T     EXCHANGE 
)-X-(IQN(l) -IQN(2)  ) 
)-X*(IQN{3)-IQN(4) ) 
=MCNU-Q\'L  ) 
E    HEIGHT 

)-X^MIWH(  1)-IWH(2)  ) 
)-X^-  (IWH(3)-IWh(4)  ) 
"( WHU-WHL  ) 
LL    HEIGHT 

)-X^(ISH{l )-ISH(2)  ) 
)-X^( ISH(3)-ISH(4) ) 
-(SHU-SHL) 
SCALING    FACTOR 


)+160. ) /5. 
SC)+160.)/5 

2. 
2. 
I  t4 
)/16. 
)/16. 
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CON 


MEC 
12 


13 

CCK 

10 

GRI 


5 

30 
31 
32 
33 
34 
35 


9999 


CP  = 

pys  = 

DCO 
IF( 
VECTIV 
DT= 
CN  = 
DCG 
H/5NIC/i 
CO  CH  = 
DM  = 
IF( 
DL  = 
GQ 
00  CL  = 
V/DI  V 
OC  A  =  U 
C  SP/^C 
G  = 
8  =  ( 
DC  = 
PPL 
WRI 
WRI 
^RI 
WRI 
WRI 
WRI 
WRI 
FCR 
FOR 
FCR 
FCR 
FCR 
FCH 
FOR 
GO 
STC 
END 


0.935 

I.  024 

N=0.  0 

GN.GE.0.0)  GO  TO  1200 

E  MIXING 

(  lOO.'^(T-TS)  )/(PLI-600.  ) 

(2  .•'■<GN)/{CP^RW^DT) 

N=(PLI'^*2  +  CM)  --^^0.5 

L  MIXING 

{WH--^''^2  +  SH=^^2)'^'^0.5 

12.5'^CH 

CM.GE.DCON)  GO  TO  1300 

DCON 

TO  ICOO 

DM 

CORRECTION 

(1)  +  U(3)-U(2)-U(4)  +  V(3)  +  V{4)-V{1)-V(2) 

ING    'G' 


02 

CC 
00 
OC 
OC 
00 
CC 


A'^0. 
B^'DL 
0  =  DL 
TE(  b 
TE(6 
TE(6 
TE{6 
TE(6 
TE(6 
TE(6 
MAT  { 
MAT( 
MAT  ( 
MAT( 
MA  T  { 
MAT  ( 
MAT( 
TO  9 
P 


864)/l4.=^G) 


+  DD 
50 
30 
31 
32 
33 
34 
35 


2)  KT,KD 
00)  DCGN 

DM 

DL 

DD 

PRLD 

RLD 
13X, 


COO 


CO) 

00) 

00  ) 

GO) 

00) 

,//// 

t23X, 

i23X, 

T23Xt 

,23Xr 

»23X, 

,23X, 


'FORECAST  FOR 
•CONVECTIVE  MLD 
•FORCED  MLD  IS  • 
' LAYER  DEPTH  IS 
•CONV/OIV  CORR  I 
•FORECAST  MLD  IS 
•ACTUAL  MLD  IS  • 


S  12,  *l 

•  ,F8.2  ) 
.2) 

8.2) 
,  F  8 .  2  ) 
•,F8.2) 
,F8.2) 


IS 
.Ffi 

■  »F 
<;  I 


12,'  JAN.«) 
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APPENDIX  I 
FORECAST    MLD    FOR    STATION    NOVEMBER 

FORECAST    FOR      OZ    13    SEP. 

CONVECTIVE    MLD    IS  0.0 

FCRCED    MLC    IS  53.71 

LAYER    DEPTH    IS  53 .71 

CONV/DIV    CGRR    IS  -0.37 

FORECAST    MLD    IS  53.34 

ACTLAL    MLD    IS         110.00 


FORECAST    FOR     12Z    13    SEP . 

CCNVECTIVE    MLD    IS  0.0 

FORCED    MLD    IS  59.25 

LAYER    DEPTH     IS  59.25 

COKV/DIV    CCRR    IS  -0.11 

FORECAST    MLD    IS  59.14 

ACTUAL    MLD    IS  60.00 


FORECAST    FCR      02    14    SEP. 

CONVECTIVE    MLD     IS  110.26 

FORCED    MLC    IS  55.67 

LAYER    DEPTH    IS  110.26 

CONV/DIV    CORR    IS  -0.17 

FORECAST    MLD    IS  110.09 

ACTLAL    MLD    IS         180.00 


FORECAST     FOR     12Z    14    SEP. 

CONVECTIVE    MLD    IS  60.21 

FORCED    MLD    I S  55. 18 

LAYER    DEPTH     IS  60.21 

CONV/DIV    CORR    IS  -0.11 

FORECAST    MLD    IS  60.10 

ACTUAL    MLD    IS  160,00 


FORECAST    FCR    12Z    15    SEP. 

CONVECTIVE    MLD    IS  160.49 

FORCED    MLC    IS  73.91 

LAYER    DEPTH     IS  160.49 

CONV/DIV    CORR     IS  1.09 

FORECAST    ^'LD    IS  161.58 

ACTUAL    MLD    I S  130.00 


FORECAST    MLD    FOR    STATION    NOV  EMBER (CONT ) 


FORECAST    FOR       OZ    16    SEP. 

CONVECTIVE    MLD    IS  150.24 

FORCED    MLC    IS  51.32 

LAYER    DEPTH     IS  150.24 

CONV/DIV    CORR    IS  1.11 

FORECAST    NLD    IS  151.35 

/iCTUAL    MLD    IS  150.00 


FORECAST     FOR    122    16    SEP. 

CONVECTIVE    MLD     IS  130.50 

FORCED    MLD    IS  53.70 

LAYER    DEPTH     IS  130  .50 

CONV/DIV    CORR    IS  1.01 

FORECAST    MLD    IS  131.51 

ACTUAL    MLD    IS         18  0.00 


FORECAST    FOR       OZ    19    SEP. 

CONVECTIVE    MLD    IS  190.11 

FORCED    MLC    IS  37.08 

LAYER    DEPTH    IS  190.11 

CONV/DIV    CORR    IS  1.03 

FORECAST    MLD     IS  191.14 

/SCTUAL    MLD    IS  20.00 


FOREC/iST    FCR    12Z    19    SEP. 

CONVECTIVE    MLD    IS  160.08 

FORCED    MLD    IS  32.67 

LAYER    DEPTH    IS  160  .08 

..   CONV/DIV    CORR    IS  0.  11 

FORECAST     MLD    IS  160.19 

ACTUAL    MLD    IS         14  0.00 


FORECAST    FOR       OZ    20    SEP. 

CCNVECTIVE    MLD    IS  22.52 

FORCED    MLC    IS  27.62 

LAYER    DEPTH    IS  27.62 

CONV/DIV     CORR    IS  0.21 

FORECAST    MLD    IS  27.83 

/ACTUAL    MLD    IS  170.00 
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FORECAST    MLC    FOR    STATICN    NOVEMBE R ( CCNT ) 


FORECy^ST     FOR     12Z    2C    SEP 


CONVECTIVE    MLD 
FORCED    MLD    I  S 
LAYER  DEPTH  IS 
CCNV/DIV    CORR    IS 
FORECAST    f^LD    IS 
ACTUAL    MLD    IS 


IS  140.22 

29.86 
140.22 

0.50 
140.72 
170.00 


**  *  :ic^  A  ^  :^  Jt  ^  ^  jt:  ^  :^  *  ~  ^  A  ^ :^  :,5t  .-X  jje  :^  :<t -; *  ^  :^ A  sjc  ^  >:x  :<<  ;?  3}:  :»£  :^  j^  :^  :^  A  :^  ^^ 
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FORECAST    MLD    FOR    STATI  CiN    NOVEMBER 


FORECAST     FOR      OZ    12    JAN. 

CGNVECTIVE    MLD    IS  550.06 

FORCED    MLD    I S  132. 37 

LAYER    DEPTH     IS  550.06 

CONV/DI V    CORR    IS  1  .90 

FORECAST    MLD    IS  551.96 

ACTLAL    MLD    IS         530.00 


FORECAST    FOR    12Z    12    JAN. 

CONVECTIVE    MLD    IS  560.04 

FORCED    MLC    IS  125.29 

LAYER    DEPTH     IS  560. 04 

CONV/DIV     CORR    IS  1.44 

FORECAST    NLD    IS         561.48 
/ACTUAL    MLD    IS         570.00 


FORECAST     FCR      OZ    13    JAN 


CGNVECTIVE    MLD 
FORCED    ML  D    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    1 
FORECAST    MLD    IS 
ACTUAL    MLD    IS 


IS  530.04 

138.08 
530  .04 
S  0.76 

530.80 
510 .00 


FORECAST    FOR     12Z    13    JAN. 

CGNVECTIVE    MLD    IS  570.02 

FORCED    MLC    IS  73.65 

LAYER    DEPTH     IS  570.02 

••  CCNV/DIV    CORR    IS  1.94 

FORECAST    MLD    IS         571.95 
ACTUAL    MLD    IS         530. CO 
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FCRECAST    MLD    FOR    STATION    PAPA 


FORECAST    FOR       OZ    L3    JAN. 

CONVECTIVE    MLD    IS  390.15 

FORCED    MLC    IS  164.92 

LAYER    DEPTH    IS         390. 15 
CONV/DIV    CQRR    IS  10.14 

FORECAST    NLD    IS         400.29 
ACTUAL    MLD    IS         360. CO 


FORECAST    FOR    12Z    13    JAN 


CONVECTIVE    MLD 
FORCED   MLD    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    I 
FORECAST    MLD    IS 
ACTLAL    MLD    IS 


IS  0.0 

130.68 
130.68 
S  -0.88 

129.80 
430.00 


FORECAST    FOR       OZ    19    JAN. 

CONVECTIVE    MLD    IS  430.82 

FORCED    MLC    IS         139.04 
LAYER    DEPTH     IS         430.82 
CONV/DIV    CORR    IS  8.90 

FORECAST    MLD     IS         439.72 
/ACTUAL    MLD    IS         390.00 


FORECAST     FCR    12 Z    19    JAN 


CONVECTIVE    MLD 
FORCED    MLD    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    I 
FCRECAST    MLD     IS 
ACTLAL    MLD    IS 


IS  382.14 

171.  10 
382.14 
S  3.53 

385.67 
410.00 


FORECAST     FOR       OZ    20    JAN. 

CCNVECTIVE    MLD    IS  391.68 

FORCED    MLC    IS         172.31 
LAYER    DEPTH     IS  391.68 

CONV/DIV    CORR    IS  10.94 

FORECAST    MLD    IS         402. 63 
/SCTUAL    MLD    IS  380. CO 
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FORECAST    MLO    FOR    STATION    PAPA(CCNT) 


FORECAST    FOR    L2Z    20    JAN. 

CCNVECTIVE    MLD    IS  410.00 

FORCED    MLD    IS  134.28 

LAYER    DEPTH     IS  410.00 

CONV/DIV    CORR    IS  9.59 

FORECAST    MLD    IS  419.59 

ACTUAL    MLD    IS         460.00 


***5»j5>«  J!^^fe*5}cA5^:5if^:J;**:^«^A*^V^yc^|S5jt5!«*^^isyt>!<::*A5h*:^:**A*>>i:^5^A  *  3i<  :<t  J;  *  3;t *  * >:<  *  i^  A  *  *  *  a  jje  *  5^  :(;  >^  ?!<  •-■:«  >!:  A  * « s^  * ' 

»  «  :fr  sit  *  i(  sli  :jt  ^e  sjt  ;it  sjs  ^  ^ -!  A  *  *  :{{ :^  :{f -i  ^  :i}[  ,^  A  ^  ^  *  rs  :p  5{:  *  ;^  *  ^!!  ^  3{c  •-!«  A  .1«  :?^ 
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FCRECAST    MLD    FOR    STATION    ECHO 


FORECAST    FOR       OZ    12    JAN 


CONVECTIVE    MLD    IS 
FORCED    MLD    IS 
LAYER    DEPTH     IS 
COKV/OIV    CORR    IS 
FORECAST    KLD    IS 
/iCTUAL    MLD    IS 


898 
83.52 
898.  16 
11.99 
910.15 
110.00 


16 


FORECAST    FOR    12Z    12    JAN 


CONVECTIVE    MLD 
FORCED    MLD    IS 
LAYER    DEPTH    IS 
CONV/DIV    CORR    IS 
FCRECAST     MLD     IS 
ACTUAL    MLD    IS 


IS  6A9.67 

92.04 
649. o7 
11.  13 
660.81 
580.00 


FORECAST     FOR       OZ     15    JAN. 

CCNVECTIVE    MLD     IS  769.89 

FORCED    MLD    IS  33.51 

LAYER    DEPTH    IS  769.89 

COKV/DIV    CORR    IS  -0.14 

FORECAST    MLD    IS  769.75 

ACTUAL    MLD    IS         620.00 


FORECAST    FCR    12Z    15    JAN. 

CONVECTIVE    MLD    IS  659.92 

FORCED    MLD    IS  52.49 

LAYER    DEPTH    IS         659.92 
CONV/DIV    CORR     IS  18.80 

FCRECAST    MLD    IS  678.72 

ACTUAL    MLD    IS         650.00 


FORECAST     FOR       OZ    16    JAN. 

CONVECTIVE    MLD    IS  0.0 

FORCED    MLD    IS  45.71 

LAYER    DEPTH    IS  45.71 

CONV/DIV    CORP    IS  1.30 

FORECAST    MLD    IS  47.01 

ACTUAL    MLD    IS  100.00 
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FORECAST    MLD    FOR    STATION    ECHGCCCNT) 


FORECAST     FOR     L2Z    16    JAN. 

CCNVECTIVE    MLD    IS  649.89 

FORCED    MLD    IS  6  8. 38 

LAYER    DEPTH     IS  649.89 

CGNV/DIV    CORR    IS  It3.86 

FORECAST    MLD    IS  665.75 

ACTUAL    MLD    IS         670.00 
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CONTINUCUS    MLD    FORECAST     FOR    STATION    NOVEMBER 


FORECAST    FOR    12Z    IL    SEP. 

CONVECTIVE    MLD    IS  0.0 

FORCED    MLD    IS  39.83 

LAYER    DEPTH    IS  39.83 

CONV/DIV    CORR     IS  -0.73 

FORECAST    MLD    IS  39.10 

ACTUAL    MLD    IS  130.00 


FORECAST     FOR     121     12    SEP. 

CONVECTIVE    MLD    IS  0.0 

FORCED    MLD    I S  49.  13 

LAYER    DEPTH     IS  49.  13 

CPiW/DIV    CORR    IS  -0.41 

FORECAST    MLD    IS  48.72 

ACTUAL    MLD     IS  160.00 


FORECAST     FCR    12Z    13    SEP. 


CONVECTIVE    MLD 
FORCEC    MLD     IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR     IS 
FORECAST     MLD    IS 
ACTUAL    MLD    I S 


IS  0.0 

59.25 
59.25 
-0.  11 
59.14 
60.  00 


FORECAST     FOR     12Z     14    SEP. 

CONVECTIVE    MLD    IS  59.36 

FORCED    MLD    IS  55.  18 

LAY  ER    DEPTH     IS  59.36 

CONV/DIV    CORR    IS  -0.11 

FORECAST    MLD    IS  59.2  5 

ACTUAL    MLD    IS  160.00 


FORECAST     FCR    122    15    SEP. 

CONVECTIVE    MLD    IS  60.86 

FGRCEl    MLD    IS  73.91 

LAYER    DEPTH    IS  73.91 

CONV/DIV    CORR    IS  0.50 

FORECAST    MLD    IS  74.41 

ACTUAL    MLD    IS  130.00 


L 
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FCRECAST    MLD    FOR    STATION    NOV EMBE R ( CONT ) 


FORECAST    FOR    12Z    16    SEP. 

CONVECTIVE    MLD    IS  75,38 

FORCED    MLD    IS  53 .70 

LAYER    DEPTH    IS  75.33 

CONV/DIV    CORR    IS  0.58 

FCRECAST    f^LD    IS  75.97 

ACTU/iL    MLD    IS  180.00 
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CONTINUCLS    MLD    FCRECAST     FCR    STATION    NOVEMBER 


FORECAST    FOR       OZ    11    SEP. 

CGNVECTIVE    MLD    IS  0.0 

FORCED    MLC    IS  29.02 

LAYER    DEPTH    IS  29.02 

COKV/DIV    CORR    IS  -0.44 

FORECAST    KLO     IS  28.57 

ACTUAL    MLD    IS  160.00 


FCRECAST     FCR      OZ    12    SEP. 

CONVECTIVE    MLO    IS  0.0 

FORCED    MLD    IS  18.86 

LAYER    DEPTH    IS  18  .86 

CONV/OIV    CORR    IS  -0.20 

FORECAST     MLD    IS  18.66 

ACTUAL    MLD    IS  120.00 


FORECAST     FOR       OZ     13    SEP. 

CONVECTIVE    MLD    IS  0.0 

FORCED    MLD    IS  53.71 

LAYER    DEPTH     IS  53.71 

CONV/DIV    CORR    IS  -0.37 

FORECAST    MLD    IS  53.34 

ACTUAL    MLD     IS  110.00 


FCRECAST     FCR       OZ    14    SEP 


CONVECTIVE    MLD 
FORCED    MLD    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    IS 
FORECAST    MLD    IS 
ACTUAL    MLD    IS 


IS 


53 
55.67 
55.67 

-C.  09 
5  5.59 
180.00 
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CONTINUOUS    MLD    FORECAST    FOR    STATION    NCVEI^BER 


FORECAST     FOR       OZ    12    JAN. 

CCNVECTIVE    MLD    IS  550.06 

FORCED    MLD    IS  132. 37 

LAYEf;    DEPTH     IS  550.06 

CCNV/DIV    CORP    IS  L.90 

FORECAST    MLD    IS  551.96 

ACTUAL    MLD    IS         530.00 


FORECAST    FOR       OZ    13    JAN. 


CONVECTIVE    MLD 
FGRCLO    MLC    IS 
LAYER    DEPTH     IS 
CO^V/C^^/    CORR    IS 
FCRECAST    PLD    IS 
ACTUAL    MLD     IS 


IS  551.99 

138  .08 
551.99 

0.79 
552.78 
510. CC 


FCRECAST     FCR      OZ    14    JAN 


CONVECTIVE    MLD 
FORCED    ML  D    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    I S 
FORECAST    MLD    IS 
ACTUAL    MLD    IS 


IS 


552 
83.  18 
552.81 

1.61 
554.42 
520.00 


61 


FORECAST    FOR       OZ    15    JAN. 

CCNVECTIVE    MLD    IS  0.0 

FORCED    WLD    IS  42.45 

LAYER    DEPTH     IS  42.45 

CCNV/DIV    CORR    IS  0.33 

FORECAST    MLD    IS  42.78 

ACTUAL    MLD    IS         540. CO 


FORECAST     FCR       OZ    16    JAN 


CONVECTIVE    MLD 
FORCED    MLD     IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    IS 
FCRECAST     MLD    IS 
ACTUAL    MLD    IS 


IS  45.77 

82.02 
82.02 
-0.  23 
8  1.79 
440.00 
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FORECAST    MLD    FOR    STATION    NO V EMSER (CONT ) 


FORECAST     FOB      02    17    JAN. 

CONVECTIVE    MLD    IS  84.  14 

FCRCEC    MLC    IS  114.14 

LAYER    DEPTH     IS  114.14 

CONV/DIV    CORR     IS  C.44 

FORECAST    f^LD    IS  114.58 

ACTLAL    MLD    IS         540. GO 


FORECAST    FOR       OZ    18    JAN. 

CONVECTIVE    MLD    IS  0.0 

FORCED    MLD    IS  46.32 

LAYER    DEPTH     IS  46.32 

COKV/DIV    CORR    IS  -0.57 

FORECAST    MLD    IS  45.75 

ACTUAL    MLD    IS         520.00 


FORECAST    FCR      02    19    JAN. 

CONVECTIVE    MLD    IS  46.  17 

FORCED    MLC    IS  77  .63 

LAYER    DEPTH     IS  77. 63 

CGNV/DIV    CORR    IS  -0.29 

FORECAST     NLO    IS  77.35 

ACTUAL    MLD    IS         500. CO 


FORECAST     FOR      0 2    20    JAN 


CONVECTIVE    I^LD    I 
FORCED    MLC    I  S 
LAYER    DEPTH     IS 
CGNV/DIV    CORR    IS 
FORECAST    MLD    IS 
/5CTLAL    MLC    IS 


77.62 
^7.75 
87  .75 
-0.26 
87.48 
590.00 
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CCMINUGUS    MD    FORECAST    FOR    STATION    PAPA 


FORECAST    FOR       OZ    L3    JAN. 

CONVECTIVE    MLD    IS  390.15 

FORCED    MLD    IS         164.92 
LAYER    DEPTH    IS  390.  15 

CQNV/OIV    CORR    IS  10.1.4 

FORECAST    NLD    IS         400.29 
ACTUAL    MLD    IS         360. CO 


FORECAST    FOR      OZ    14    JAN 


CONVECTIVE    MLD 
FORCED    MLD    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    I 
FORECAST    MLD    IS 


IS  0.0 

231.  18 
231.18 
S  1.57 

232.75 


ACTUAL    MLD    IS      9999.00 


FORECAST    FOR       OZ    15    JAN. 

CONVECTIVE    MLD    IS  0.0 

FORCED    MLD    IS         332.13 
LAYER    DEPTH    IS         332.  13 
CCNV/DIV    CORR    IS  4.07 

FORECAST    NLD    IS         336.20 
ACTUAL    MLD    IS       9999. CO 


FORECAST     FCR      OZ    16    JAN 


CONVECTIVE    MLD 
FORCED    MLD    IS 
LAYER    DEPTH     IS 
CONV/DIV    CORR    I 
FORECAST    MLD    IS 
ACTUAL    MLD    IS 


IS         337.39 
211.43 
337.39 
S  -1.25 

336.  14 
330.00 


FORECAST    FOR       OZ    17    JAN. 

CONVECTIVE    MLD    IS  338.42 

FCRCEO    MLD    IS         204.18 
LAYER    DEPTH    IS  338.42 

CONV/DIV    CORR    IS  1.04 

FORECAST    MLD    IS         339.46 
ACTUAL    MLD    IS       9999.00 
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FORECAST    NLD    FOP    STATION    PAPA{CONTJ 


FORECAST    FOR       OZ     18    JAN. 

CCNVF.CTIVE    MLD    IS  341.38 

FORCED    MLC    IS         192. C8 
LAYER    DEPTH     IS  341.38 

CGNV/DIV    CORR    IS  6.54 

FORECAST    NLD    IS         347.92 
ACTUAL    MLD    IS         430.00 


FORECAST    FCR      OZ    19    JAN 


CONVECTIVE    MLD 
FORCED    MLC    IS 
LAYER    DEPTH    IS 
CONV/DIV    CORR    I 
FORECAST    MLD    IS 
ACTLAL    MLC    IS 


IS  349.43 

L39,04 
349.43 
S  7.22 

356.65 
390.00 


FORECAST     FOR       OZ    20    JAN 


CONVECTIVE    MLD    IS  358.7 

FORCED    MLD    IS  172.31 

LAYER    DEPTH     IS  358.78 

CGNV/DIV    CORR    IS  10.02 

FORECAST    MLD    IS  368.80 

ACTUAL    MLD    IS  380.00 


8 


1 
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